UNIVERSITY of 
NOTRE DAME 


NASA/USRA UNIVERSITY 
ADVANCED DESIGN PROGRAM 
1989-1990 

UNIVERSITY SPONSOR 
BOEING COMMERCIAL AIRPLANE COMPANY 


FINAL DESIGN PROPOSAL 


THE DRAG-n-FLY 

A Proposal in Response to a Low Reynolds Number 
Station Keeping Mission 

May 1990 


Department of Aerospace and Mechanical Engineering 
University of Notre Dame 
Notre Dame, IN 46556 


(NASA-C«-18<.b6o) LRAG-n-FLY: A PROPOSAL IN NPG-25130 

RESPONSE TO A OOw REYNOLDS NUNbER STATION 
KEEPING MISSION Pin?! Design Proposal 

(Notr^ o d Tit: Univ.) i- 7 3 p CSCL OIC unci as 

03/05 0769191 


‘Drag - n - ^Cy 


A Proposal in Response to a Low Reynolds Number Station Keeping Mission 




John Nieliaus Joey Brenner 

Jenny Neumann Peter Gendron 

Pat Deviny 


Submitted May 7, 1990 


Table of Contents 


EXECUTIVE SUMMARY 

1 PRELIMINARY CONCEPT DESIGN 

1 1 MISSION REQUIREMENTS AND OBJECTIVES 

1.1.1 MISSION DEFINITION 

1.1.2 AIRCRAFT DESIGN 

1.1.3 DESIGN OBJECTIVES 

1.1. 3.1 TAKE-OFF 

1. 1.3.2 CRUISE 

1. 1.3.3 TURNS 

1. 1.3.4 FLIGHT PATH /LANDING 

1. 1.3.5 GENERAL CONFIGURATION REQUIREMENTS 

1.2 CONCEPT SELECTION 

2 CONCEPT DESIGN 

2.1 AERODYNAMIC DESIGN 

2.2.1 WING DESIGN 

2.1.2 AIRFOIL SELECTION 

2. 1.2.1 AIRFOIL THICKNESS 

2. 1.2.2 AIRFOIL SELECTION DATA AND CRITERIA 

2.1. 2.2.1 STALL BEHAVIOR 

2.1. 2.2.2 STALL ANGLES 

2.1. 2.2.3 AIRFOIL SHAPE 

2.1. 2.3 THE SPICA AIRFOIL 

2.1.3 DRAG ESTIMATION 

2.1. 3.1 PARASITE DRAG 

2. 1.3.2 INDUCED DRAG 

2.2 PROPULSION 

2.2.1 SYSTEM SELECTION 

2.2.2 PROPELLER SELECTION 

2.2.3 ENGINE CONTROL 

2.3 STABILITY AND CONTROL 

2.3.1 STATIC MARGIN 

2.3.2 LONGITUDINAL STABILITY 

2.3.2. 1 HORIZONTAL TAIL SIZING, LOCATION, AND 
ORIENTATION 

2.3.2.2 FINAL DESIGN 

2.3.3 ROLL AND DIRECTIONAL STABILITY AND CONTROL 

2.3.3. 1 ANALYSIS TOOLS 

2.3.3.2 OBJECTIVES AND CONSTRAINTS 

2.3.3.3 DIHEDRAL, RUDDER AND VERTICAL TAIL SIZING 

2.4 STRUCTURAL ANALYSIS 

2.4.1 V-N DIAGRAM 

2.4.2 STRUCTURAL COMPONENTS 

2.4.2. 1 WING 

2.4.2.2 FUSELAGE 

2.4.2.3 EMPENNAGE 

2.4.2.4 LANDING GEAR 

2.4.2.5 MATERIALS SELECTION 

2.4.2.6 GENERAL CONSTRUCTION 

2.4.2.6.1 INTERNAL LAYOUT 

2.4.2.6.2 SCALE RPV DRAWING 

2.5 WEIGHT AND CENTER OF GRAVITY ESTIMATION 


1-1 

1-1 

1-1 

1-1 

1-1 

1-2 

1-2 

1-3 

1-4 

1-4 

1-4 

2 . 1-1 

2.1-1 

2.1-1 

2.1-4 

2.1-4 

2.1- 4 

2.1- 5 

2.1- 5 

2.1- 7 

2.1- 7 

2 . 1 - 8 

2.1- 9 
2 . 1-10 

2 . 2 - 1 
2.2-1 

2.2- 3 

2.2- 4 

2.3- 1 

2.3- 1 

2.3- 3 
2.3.4 

2.3- 8 

2.3- 9 

2.3- 9 

2.3- 10 

2.3- 11 

2.4- 1 
2.4-1 
2.4-1 

2.4- 1 

2.4- 8 

2.4- 8 

2.4- 11 

2.4- 17 

2.4- 18 

2.4- 18 

2.4- 19 

2.5- 1 



2.5.1 COMPONENT WEIGHTS 

2.5.2 CENTER OF GRAVITY ESTIMATION 

3 PERFORMANCE ESTIMATION 

3.1 TAKE-OFF AND LANDING ESTIMATES 

3.2 RANGE AND ENDURANCE 

3 3 POWER REQUIRED AND POWER AVAILABLE 
3.4 CLIMB AND GLIDE PERFORMANCE 

4 TECHNOLOGY DEMONSTRATOR 

4.1 FINAL CONFIGURATION 

4 2 TECHNICAL DEMONSTRATOR FLIGHT TEST 

5 HIGH ALTITUDE FEASIBILITY 


2.5- 1 

2.5- 3 
3-1 
3-1 
3-2 
3-3 

3- 5 

4- 1 
4-1 

4- 2 

5- 1 



EXECUTIVE SUMMARY 


The Drag-n-Fly is a remotely piloted, low Reynolds number vehicle. It was designed to 
maintain level controlled flight and fly a closed course at flight speeds corresponding to 
Reynolds numbers of less than 2 x 10 5 and as close to 1 x 10 5 as possible. The success of 
the mission will be associated with achieving the lowest mean chord Reynolds number 
possible and maximizing loiter time cm the course. 

A three-view drawing and specifications summary follow this executive summary. 

The flight plan for the Drag-n-Fly calls for the vehicle to ascent to a cruise altitude of 25 ft. 
Once achieved, the Drag-n-Fly will demonstrate it’s maneuverability by flying within a 
restricted altitude range on a figure "8” course with a spacing of 150 feet between the 
markers which define the course. The Drag-n-Fly will complete three laps and then a final 
oval to bring the RPV back around in preparation for landing. 

The airfoil selected for the Drag-n-Fly is a Spica chosen for its high lift coefficient at low 
Reynolds number. The actual lifting surface for the Drag-n-Fly is a 8.5 foot long wing 
with total surface area of six square feet and aspect ratio of 12. There is no sweep or twist 
associated with the wing and the taper ratio is 1 .0. During the flight the lifting surface 
expects a wing loading of approximately 7.1 Oz/ft 2 . 

The propulsion system for the Drag-n-Fly consists of a 10" diameter propeller mounted on 
the front of the vehicle. The propeller, the Zinger 10-6, is driven by the ASTRO 05 electric 
motor and eight 500 MAH Nickel Cadmium batteries. This motor/battery combination was 
selected not only because it is capable of providing the thrust needed to accomplish the 
mission, but also because of its lightweight nature. An electronic speed control will also be 
utilized to maintain altitude through the turns. 

Structural support for the Drag-n-Fly comes from four box beams running the length of the 
fuselage. These box beams will be constructed of Balsa wood. The strongest part of the 
aircraft is the front fuselage, since the motor and avionics are located in this region. This 
area will be reinforced by panels to increase the strength of the front fuselage. The spars of 
the wing will be made of spruce and the ribs will be formed of Balsa wood. 



The tail and horizontal stabilizers are located for aft of the lifting surface in order to assure 
proper static stability. The vertical tail area is 0.5 square feet and the horizontal tail area is 
1.05 square feet. Two movable control surfaces are used to maneuver the Drag-n-Fly 
during the figure eight course. A rudder will be implemented to control yaw and an 
elevator to control the pitch during the flight course, and both control surfaces will be 
actuated by mirco-servos. The roll stability will be handled by dihedral angle of 10.5 
degrees. 

The present design for the Drag-n-Fly will meet the criteria for the present mission. Some 
areas of concern are in the construction of the dihedral angle on the wing when it is 
mounted to the fuselage, control of the aircraft once in flight (will the control surfaces 
deflect enough to maneuver the aircraft), and the fact that no one in Design Group A has 
any prior experience in building. While there are areas of concern, this design group is 
very confident that this aircraft will be a highly-maneuverable RPV with the capability to fly 
at low Reynolds number regimes. 
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1 preliminary concept design 


1.1 MISSION REQUIREMENTS AND OBJECTIVES 

1.1.1 MISSION DEFINITION 

The objective of this mission is to develop a remotely piloted vehicle (RPV) which will 
maintain level controlled flight and fly at speeds corresponding to Reynold's numbers less 
than 2 x 10 5 and as close to 1 x 10 5 as possible. This mission is defined with the following 
requirements and constraints. The RPV is required to fly a closed figure eight pattern 
within a restricted air space. The course envelope is constrained by the dimensions of the 
flight arena (Loftus Center) and two pylons spaced 150 feet apart. The maximum 
allowable altitude is 25 feet, and the RPV must never exceed this height at any point on the 
course. The RPV must take-off from the pit area, complete three figure eight laps, and 
finally, return to the pit area for landing. The greatest measure of merit is awarded for 
accomplishing this mission with the lowest mean chord Reynold’s number in conjunction 
with maximizing the loiter time on the closed course. 




Maintain level flight at 1 x 10^ < Re < 2 x 10-*. 
Fly within a restricted airspace, (Loftus Center). 
Follow a closed figure eight pattern. 

Never exceed altitude of 25 feet 
Take-off and land in pit area. 


1.1.2 AIRCRAFT DESIGN 

The proposed aircraft design must satisfy the following requirements and constraints. Only 
those propulsion systems, which have non-airbreathing engines and do not emit any mass, 
may be used. The RPV must take-off from the ground, and land on the ground, and do so 
within the designated 150 foot pit area. It must be able to be remotely controlled, with 
ample maneuverability and control, so that it can be flown within the limited airspace by a 
pilot with minimal flying experience. The aircraft must be of adequate strength and size in 
order to carry An instrument package payload which has a weight of two ounces and 
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dimensions of 2” x 2" x 2". In addition, the complete aircraft must be able to be 
disassembled to fit within a 2' x 2’ x 4’ container for transportation and storage. 




• Non- Airbreathing engines only. 

• Take-off and land from the ground. 

• Must be remotely controlled. 

• Fit in a 2’x 2'x 4' storage container when dismantled. 

« Carry a 2 ounce , 2"x 2"x 2" payload. 


1.1.3 DESIGN OJECTIVES 

The mission and system were further fixed by the design objectives and goals set by the 
design team. These objectives and goals were determined independent of the RPV 
configuration. 

1.1.3.1 TAKE-OFF 

The first phase of the mission is the take-off and climb stage. The main objective of the 
take-off phase is to have the airplane achieve lift off and reach a cruise altitude of 15 feet by 
the end of the 150 foot pit area. This goal, in turn, determined the ground roll distance and 
ascent angle. In order to reach the desired altitude of 15 feet by the end of the 150 foot pit 
area, a maximum ground roll distance of approximately 30 feet with an ascent angle of 10 
degrees and a rate of climb of approximately 3. 1 feet per second is necessary. The landing 
gear will be a permanently attached system with no assisting launch device being 
implemented. Wheel control will not be used since the rudder will provide sufficient 
control to steer the plane whife taxiing. 

1. 1.3.2 CRUISE 

Once altitude has been reached, a banking maneuver must be performed immediately, in 
order to perform the required turn. A large turn radius is desired to maximize the loiter 
time through a turn. This corresponds to a longer flight distance as well as a lower bank 
angle. The lower bank angle reduces the power needed to propel the RPV through each 
turn. The targeted radius is 55 feet, allowing five feet on the outside to compensate for 
pilot error and response time. The turn will encompass 190 degrees. The remaining 
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portion of the flight lap consists of two straight flight 

the two turns as shown in figure 1.1-1. 

FLIGHT TEST ARENA 


sections of equal length, connecting 
LEGEND 



40 yds 


q - Pylon 

= - Takeoff/pit area 

— - Flight path 

— - Take-off path 

— - landing path 


FLIGHT TEST COURSE 



Figure 1.1-1 


1. 1.3.3 TURNS 

During each turn, the RPV will bank at an angle of approximately 20 degrees and 
experience a wing loading of 7.1 ounces per square foot. It is desired to roll the RPV to 
the necessary bank angle in a minimum of two seconds, requiring a roll rate of 10 degrees 
per second. Turn maneuvers will be performed by deflecting the rudder and elevators. 
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Two control channels will be used to control these deflections. A third control channel will 
employ a speed control device to ensure level turns throughout the flight, as a result, the 
speed will not remain fixed during the flight. Although the speed will vary, a target value 
of 25 feet per second was set as the cruise velocity. 

1. 1.3.4 FLIGHT PATH / LANDING 

The RPV flight path for the entire mission includes a 30 foot ground roll, a 10 degree 
ascent to an altitude of 15 feet, three figure eight laps, and a final sweep around the course 
to bring the RPV into position for landing. The total distance of this flight corresponds to a 
range of 4830 feet and an endurance of 193 seconds. The RPV will descend as it 
approaches the pit area for landing. The descent angle will not exceed 10 degrees, or go 
below the minimum glide angle of three degrees. Upon touchdown, the engine will be shut 
off. It will be possible to switch the engine on and off via a remote control. Since no 
brakes will be used, a landing distance of 64 feet is needed to bring the plane to a complete 

stop. 

1. 1.3.5 GENERAL CONFIGURATION REQUIREMENTS 

The final area of design to be addressed are general configuration requirements which will 
affect any particular design chosen. The RPV weight shall remain constant during all 
phases of the flight with the target weight being three pounds. A maximum Uft coefficient 
of 1.2 is desired. A color scheme will be chosen such that the RPV will be highly visible 
to the pilot and design team, but basic, so as not to require an additional weight penalty. 


1.2 CONCEPT SELECTION 

Initially, from the seven concepts that were produced by each member there were two basic 
types of aircraft that were being considered to complete the mission. The first was an 
unconventional tandem/joined-wing aircraft with the wings joined at a point about three- 
quarters out on the semi-span. The other six being the more conventional wing-tail- 
fuselage type aircraft with the primary differences in the tails. These seven concepts were 
reduced to three such that they could be compared and one could be selected for the design. 
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The first candidate was the tandem/joined-wing aircraft as shown in Figure 1.2-1 a,b, and 
c. The front end of the fuselage was expected to be quite large tapering down to a small 
circular back end. The propeller was placed at the rear end of the aircraft creating a pusher 
type aircraft. One possible advantage of this aircraft was a reduction in weight because the 
joined wing reduces the wing bending and because only two channels of control were 
needed. Also, the double wing would disperse the surface area necessary for flight 
allowing decreases in weight possibly leading the aircraft to fly at slower speeds which 
would provide merit for this mission. Another appealing factor of this aircraft was the 
originality and uniqueness of designing and building a different type of aircraft. 

The second concept which was considered was a conventional aircraft with a high wing, 
low horizontal tail and twin vertical tails; one each on the tips of the horizontal tail. This 
aircraft would have three channels of control; rudder, elevator and speed. Roll maneuvers 
would be achieved from a combination of the rudder and a dihedral angle. The anticipated 
advantage of the twin vertical tails was the extra rudder control power (without an 
excessively large rudder ) needed to travel the figure eight course with no ailerons. For 
ease of construction, the fuselage cross section was square in design and the wing had no 
taper. In addition the wing had a high aspect ratio for better lift and drag characteristics. 

The third concept was essentially the same as the second concept except for the tail design. 
Only one vertical tail was proposed in this design, but again, the horizontal tail would be a 
low tail. The other primary difference being that this aircraft would have ailerons and thus 
require one more channel of control. One further difference was that the fuselage was 
round in order to decrease body drag. This concept included wing taper for reducing drag. 

To help in the selection process, the two more conventional designs were combined into 
one concept boasting some advantages of each and then this concept was compared with 
the tandem wing design. The conventional aircraft that was considered is shown in Figures 
1.2-2 a,b, and c. It included a high, tapered wing with dihedral. Taper was considered to 
decrease drag and the dihedral and high wing were presumed necessary to achieve roll 
maneuvers and increase roll stability. The wing had a high aspect ratio ( about 9) for better 
lift and drag characteristics. The horizontal tail was low to keep it out of the wake of the 
wing and there were two vertical tails to provide sufficient rudder size for roll and yaw 
control and sufficient tail area for acceptable roll stability. There were three channels of 
control, speed, rudder and elevator. Roll maneuvers were expected to be achieved through 
a healthy dihedral angle and rudder deflection. The fuselage was rounded with a three inch 
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Figure 1.2-1 a 



Figure 1.2-1 b 


SIDE VIEW 
Figure 1.2-1 c 
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diameter in the forward section, just large enough to hold the payload. It tapered down to 
one inch at the tail. The fuselage size was minimized and its shape rounded in order to 
reduce drag. The proposed landing gear was tricycle type. 

Advantages of the conventional type aircraft included the following. 

• There was a large database on this type of aircraft giving typical trends for 
RPV aircraft based on reference areas. 

• With a large data base the design team could concentrate on problems 
associated with low Reynolds number and not primarily on problems dealing 
with a new type of design. 


The other concept competing with the semi-conventional RPV was the tandem/joined- wing 
plane. This concept is shown in Figures 1.2-1 a,b, and c as was noted previously. As 
shown in the diagram, the wings are connected at about three-quarters of the semi-span of 
the main wing. The secondary wing includes what would be called ailerons but they could 
act as ailerons, rudder, and elevator since they are at an angle with reference to the horizon. 
There is also a small rudder as the anticipated yaw control with the ailerons was not 
assumed sufficient. The ailerons would be one channel of control, although it must include 
a mixer, and the rudder another channel. The fuselage is large in the front since it was 
anticipated that most of the payload and batteries could be stored in this area. The landing 
gear was again tricycle type. The advantages of this concept are as follows. 

• A possible reduction in weight of the aircraft. 

• Large wing areas, including both wings, without large, heavy, and multiple spars. 
The weight reduction could not be proven but it was speculated that the joined tandem 
wings could strengthen the aircraft in the wings and back end of the fuselage and as a result 
less structural support would be needed allowing a reduction in weight. One further 
appealing factor of this concept was the originality. 

With the advantages came some disadvantages of the tandem/joined-wing design. 

• Their was no database for this type of aircraft. 

• Control surface mechanisms may be difficult to integrate into the design. 

• The tandem/joined-wing configuration may pose difficulties in longitudinal static 
stability due to the odd horizontal tail or second wing. 

The lack of a database would make it difficult to validate a lot of the results. The 
difficulties with the control mechanisms may be just one of the many construction 
difficulties in this design. 


1-9 



Measuring the advantages and the disadvantages, the semi-conventional wing, fuselage, tail 
aircraft was selected. The primary reasons being the data base availability and the 
uncertainty of whether speed and weight could be reduced with the tandem wing aircraft. 
Structural study and research of the tandem/joined-wing aircraft produced no positive 
results as the whether the advantages could actually be materialized. Also, the advantage of 
a data base for the concept selected made the initial design phase simpler allowing more 
time for the higher level design. 

Several changes in the original concept took place as more detailed analysis was done. 
Reducing the twin vertical tail to a single vertical tail was one major change. The necessary 
size of the vertical tail for roll control was not as large as anticipated. In addition, attaching 
twin vertical tails would have required adding structure outboard on the horizontal tails thus 
adding weight. Another change was the the size of the fuselage. As a truss style fuselage 
was decided upon to minimize weight, the height of the fuselage was increased. The 
increase created a larger moment of inertia, thus lowering the bending moments. The 
height changed from 3 to 6 inches and the width to 3.5 inches at the wing, tapering down 
to a 2 inch width and height at the tail. To simplify construction of the truss fuselage, the 
shape was changed from round to rectangular. Another simplification to ease construction 
was the elimination of taper of the wing. Studies on the effect of taper on the program 
LinAir demonstrated only small drag decreases that would be non-rewarding for the time 
spent making every airfoil of the wing a different size. The aspect ratio was also increased 
to even a higher value of 12. This change allowed more surface area without increasing the 
mean aerodynamic chord and Reynolds number. Despite the large increase in the span of 
the wing, the aircraft roll rate was still adequate and the structural support was still expected 
to be sufficient to hold the wing. The plane could then produce three pounds of lift with a 
coefficient of lift of 0.67. The current status of the design is as shown in Figures 1.2-3 

a,b, and c. 
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SIDE VIEW 


Figure 1.2-3 
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2 CONCEPT DESIGN 

2.1 AERODYNAMIC DESIGN 

This section discusses the wing and airfoil selection, and drag estimation method used for 
the DRAG-n-FLY. Included in this section is the lift curve for the aircraft, as well as its 

drag polar. 

2.1.1 WING DESIGN 

The wing sizing and design was one of the more flexible areas in the overall design. 
Initially, a rectangular, one-piece, high aspect ratio wing was selected for the RPV. The 
original surface area estimates came from an aircraft weight estimation, and the lift required 
to support that weight both at takeoff and at cruise. A high aspect ratio was chosen to have 
the necessary lifting surface area while minimizing the wing chord (and thus, the minimum 
chord Reynolds number). 

The selected airfoil is the SPICA, chosen from a number of other airfoils as the subject of 
an airfoil trade study. Once the airfoil was chosen, preliminary estimates of the available 
lift pushed surface area and wingspan up. The propulsion team, analyzing power required 
and power available curves, determined that a higher surface area was needed. Table 1 
demonstrates the jump in sizing over the course of the design period, from February 1 to 
March 22, 1990. 

Table 2.1-1 Changes in Wing Sizing 


Parameter 

2/1/90 

3/1/90 

3/22/90 

final concept 

AR 

9.0 

10.0 

12.0 

12.0 

S(ft2) 

3.4 

5.0 

6.0 

6.0 

b (ft) 

5.4 

7.07 

8.5 

8.5 


Several wing parameters remained unchanged from the start A taper study showed that the 
positive effects of taper, such as decreased drag, were not substantial enough to outweigh 
the difficulties in constructing a wing with rib lengths changing with spanwise distance, a 
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task requiring accurate placement and the ability to hold it all in place. Twist was never 
added for similar construction difficulties. Sweep is not needed at low speeds. 

Dihedral, however, was added to make effective turning possible without the use of 
ailerons. Rather than set the entire wing with dihedral, two outboard sections with dihedral 
were built This part of the design was driven by control analysis. 

The final wing configuration is shown in Figure 2.1-1 and has the following characteristic 
as seen in Table 2. 1-2. 

Table 2.1-2: Final Wing Characteristics 

Surface Area ,S 
Wingspan, b 
Chord, c 
Aspect Ratio, AR 
Airfoil section 
Incidence angle 
Dihedral: 

Taper 
Twist 


6 ft 2 
8.5 ft 
8.5 in 
12.0 
SPICA 
5° 

0° from fuselage to 1.2 ft along the span 
12° from 1.2 ft spanwise to tip 
None 
None 
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Figure 2.1-1 



3-D View of Left Wing 


2.1.2 AIRFOIL SELECTION 


Due to the complications associated with low Reynold's flight, such as laminar flow 
separation, selecting the proper airfoil became a crucial part of the system design. As a 
result, a study was performed to determine which airfoil section best met the needs of the 
mission. Since the RPV will fly at a low Reynold’s number (Re=100,000), only those 
airfoils that performed at well in low Re flight were considered. This, however, still left 
several airfoil sections to choose from. To solve this, a list of criteria was established in 
order to select the proper airfoil from those available. The criteria took into account the 
airfoil thickness, aerodynamic characteristics, and the shape of the airfoil. 

2.1.2. 1 AIRFOIL THICKNESS 

Thickness is an important factor in the selection of the airfoil for two reasons, strength and 
weight. Since a cantilevered wing is desired, the airfoil section needs to be thick enough 
(at least 10 % of the chord 1 ), to adequately support the main wing spar without the use of 
external bracing. However, because a low weight vehicle is desired, the airfoil section 
should not be so thick that it adds unnecessary weight. It has been found that the structural 
support gained by sections of thickness greater than 15% is very small compared to the 
additional weight gains. 2 Therefore, only those airfoils that fell in the medium thickness 
range, or those with a thickness of 10 to 14 percent of the chord, were considered. 

2.1.2.2 AIRFOIL SELECTION DATA AND CRITERIA 

Once the thickness range was selected, aerodynamic characteristics of the remaining airfoils 
were studied and compared. This was accomplished using experimental test data given on 
airfoils. The mission goal of a low Reynold's number flight, however, limited the data that 
could be used, because, as the Reynold's number decreased, so did the airfoil 
performance. At low numbers, (such as 100,000), it becomes difficult to keep the flow 
attached to the airfoil, and as a result laminar flow separation occurs. Because of this, the 
only data used, was that which was obtained at a Reynold's number similar to what was 
trying to be achieved. 


1 Korff, Walter H„ I.ifhtnlane Design. C. A. Steele Publishers, Ml Vernon, Illinois, 1935. 

2 Ibid. 
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Using lift curve slopes and drag polars 3 , the airfoils were examined for high lift (CLmax). 
the lowest drag, and a gentle stall behavior (as opposed to an abrupt stall), and a high stall 
angle (oicLmax). according to the criteria below. The design team felt that this criteria was 
the minimum the airfoil would need, based upon the mission and design objectives. In 
addition, the data was also studied for separation bubbles. Only those airfoil sections with 
the best combination of the criteria listed below, and those that did not appear to have a 
separation bubble, were examined further. 


AIRFOIL SELECTION CRITERIA 

• C L max >10 

* Low drag 

• Gentle Stall behavior 

* OCLmax> 10° 

» No separation bubbles 


2.1.2.2.1 STALL BEHAVIOR 

Although a low drag is desired, it is also important to keep stall behavior in mind. For 
instance, some of the lowest drag airfoils have the worst stall characteristics. Even though 
it may have the lowest drag, an airfoil that stalls abruptly is undesirable because, with such 
an airfoil, lift is lost without warning when the wing stalls. As a result more altitude is lost 
in trying to recover after the stall, then might have been gained from the marginally lower 
drag. 

2.1.2.2.2 STALL ANGLES 

Along with a gentle stall, airfoils with large stall angles were desired. Large stall angles 
allow for larger margins of safety which are desirable in the event that one of the 
parameters should vary from what is expected. Suppose, for instance, the friction from the 
ground causes the RPV to roll for take-off at a slower velocity than expected. The lower 
velocity would require a higher Cl value in order to achieve the lift necessary for take-off. 
For example, using the equation below, a three pound aircraft (W), with a six square feet 
wing area (S), would need a lift coefficient of about 0.87, to take-off at 20 ft/s (V). 


3 Selig, Michael S., Donovan, John F., and Fraser, David B., Airfoils at Low Speeds. H. A. Stokely, 
Publisher, Virginia, 1989. 
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L=W = 1/2 p V 2 S C L 


If the ground was so rough that it only allowed the plane to reach a velocity of 20 ft/s, the 
lift coefficient would need to be at least 1.0 to achieve lift-off. An increase in Cl, such as 
this, could put the plane at, or very near its CLm« value for some of the airfoils, as can be 
seen in Figure 2. 1 -2. Because stall occurs near Cl™*, the plane might be trying to take-off 
at its stall velocity, and as a result, it will never lift-off. In this case, however, an airfoil 
with a large stall angle would have a better chance at still being able to take off. 


2-D LIFT CURVE SLOPE FOR FINITE WINGS 
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Figure 2.1-2 

The graph above shows the lift curve slopes for the final five airfoil sections considered. 
The data shown in this plot is the for a two dimensional finite wing. It was found by 
correcting the infinite wing data for aspect ratio (AR =12) using the equation below. 4 

a = ao / [ 1 + ao/fa AR e)] 



4 Anderson, John D„ Jr., Introduction to Flight . McGraw-Hill, Inc., U. S. A. 1985. 
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The infinite data needs to be corrected because the slope of the lift curve for a finite wing 
(a), is less than that for an infinite wing (ao). The span efficiency factor (e), used in the 
above calculation was 0.87. A further discussion on how the span efficiency was 
estimated can be found in section 2.1.3. 

2.1.2.2.3 AIRFOIL SHAPE 

The final criteria used to select an airfoil was its shape. In order to make construction of 
the wing as easy as possible smooth and simple shapes were chosen. Those sections with 
complicated shapes (such as those with high camber or very thin trailing edges), were 
excluded. As a result of this evaluation, the SPICA airfoil was chosen. 

2. 1.2.3 THE SPICA AIRFOIL 

The SPICA is a flat bottom airfoil, with a thickness of 11.72 % (of the chord), and a 
camber of 4.75 % (See Figure 2.1-3). The high lift coefficients produced, combined with 
its gentle stall behavior, improve its chances for tight, low level turns. After the data was 
corrected for the aspect ratio, it was found that the SPICA has a zero lift angle of -2.0 
degrees, and could achieve a CW of 1.2 at about a 14 degree angle of attack, as shown in 

Figure 2.1-2. 



• 11.74 % thickness • CLmax =1.2 • (Xl=o- 14° 

• 4.47 % chord • cxcLmax = *2.0° 


Figure 2.1-3 

The three dimensional lift generated by the wing was produced using the lifting line 
program 5 found in Appendix C. Assuming that the lift from the fuselage and tail surfaces 
would not be significant compared to the that of the wing, the lift generated from the wing 


5 Nelson, Robert C., Atmospheric Flight Mechanics , University of Notre Dame, Indiana, 1989. 
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was taken to be the lift of the entire airplane. The lift at various attack angles was found 
using the program. This data was then used to create the three dimensional lift curve slope 
shown in Figure 2.1-4. From this graph the attack angles needed for a specific lift 
coefficient were estimated. 


3-D LIFT CURVE SLOPE FOR THE SPICA AIRFOIL 
(Lift Curve for the DRAG-n-FLY) 



Alpha [dag] 

Figure 2.1-4 

It is important to note that the stall data, (given in dashed lines) in Figures 2.1-2 and 2.1-4 
was not calculated as was the data given by the solid lines. Since the correction methods 
used (aspect ratio and lifting line program), are based upon the constant slope portion of the 
infinite wing lift curve, the stall behavior for the two and three dimensional finite wings 
could not be found using these methods. In spite of this, an estimation of stall behavior 
was made based on the data given for the infinite wing. 

2.1.3 DRAG ESTIMATION 

Overall drag is one of the key parameters limiting the design and performance of the 
aircraft. In spite of its importance the techniques used to estimate drag are approximate at 
best. The following areas are discussed: usefulness of drag analysis; estimation of parasite 
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drag; estimation of induced drag; drag polar approximation. Included in the discussion are 
the methods used in the analysis. 

In general, the primary force resisting the passage of an aircraft through the atmosphere is 
drag. The effects are felt in every aspect of design: power plant selection, wing sizing, 
and stability and control. Thus one of the primary concerns is building an aircraft that 
minimizes the negative effects of drag. Drag directly determines the power required to 
remain in the air. Basic power requirements motivate the selection of propeller and engine 
which must perform through takeoff loads, cruise, and landing. Power requirements also 
limit range and endurance. Wing sizing and airfoil selection come into play, the price of 
lift is a component of drag induced by the lift. 

2. 1.3.1 PARASITE DRAG 

The drag estimate is broken down into two areas: parasite drag and induced drag. Parasite 
drag, also known as profile drag, is the airplane drag at zero lift. It represents the effects of 
shear stress at the surface over the entire plane, the fluff body drag of a body moving 
through a fluid, boundary layer separation, and other effects. Clearly, since an analytical 
solution to the problem is unavailable (if not impossible), the parasite drag must be 
estimated empirically. One method analyzes the plane by dividing the plane into its 
components. 1 The drag on each component is then estimated experimentally. The parasite 
drag coefficient can then be estimated by summing the individual parts using 

Cdo=£ (Coi Ai / S re f) 

S re f is taken to be the wing surface area. Ai is the area on which the individual drag 
coefficient is based. Table 2.1-3 presents the drag estimated for each part. The end result 
is a value for the parasite drag coefficient, as shown below. 

Cdo=0018 

This value compares favorably with the estimates from previous year’s RPVs. 


1 Nelson. Robert C.. Atmospheric Flight Mechanics. University of Notre Dame. Indiana, 
1989. 
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Table 2.1-3: 

Parasite Drag component contributions 

Component 

cm 

Ai 

Wing 

0.007 

Wing planform area 

Fuselage 

0.110 

Max frontal area of fuselage 

Tailplane 

0.0080 

Horiz. or V re f area 

Interference 

Add 5% to Cdo 


Roughness/ 

Add 10% to Cdo 


Protuberances 



Landing gear 

Add 0.0170 



2.1.3.2 INDUCED DRAG 

The induced drag would occur even if fluids were fnctionless, since it is a component of 
the inclined lift vector. Cot, the induced drag coefficient, is approximated using the 

following equation. 


C D i = C L 2 /(*ARe) 


The induced drag is a function of the coefficient of lift, the aircraft aspect ratio, the Oswald 
efficiency factor. The induced drag varies throughout the flight, changing as the square of 
the lift coefficient. Cl is known from airfoil selection and angle of attack (as was 
discussed in the previous section 2.1.2). The only unknown parameter is e, the Oswald 
efficiency factor. This factor must also be determined from experimental data. The method 
is based on a breakdown of the contribution of the wing, the fuselage, and a miscellaneous 
term for the rest of the body. 2 The efficiency is given by the following equation. 


l/e= l/e w ing + 1/ebody + l/ e other 

The wing efficiency contribution can be found from experimental data correlating the aspect 
ratio and the wing efficiency. Using Figure 2.1-5, this value (e w ing) was found to be 
0.0.75 when the aspect ratio (AR) is equal to 12. 


2 Jensen. Daniel T., QW Predictio n Methodology for Low Reynold Nl B TlhCT Fl ig ht 
Vehicles. University of Notre Dame, Indiana, 1989. 
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1.0 



AR 

Figure 2.1-5 

The fuselage contribution can be found from the body efficiency factor, Efuseiage- It is the 
body efficiency (ebody). scaled by the cross-sectional area of the body, relative to the wing 
area. The value of Ef use iage> which equals 0.6, was found using Figure 2.1-6, with the 
aspect ratio equal to 12. 



ebody = ^-fuselage S^f / Sbody 

Once Ef US eiage is found, ebody can then be calculated using the equation above, where S re f 
is the wing area, and Sbody the maximum fuselage cross sectional area. Knowing all of the 
efficiencies, the Osward efficiency (e) was then found using the equation aforementioned. 
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e = 0.868 

From the estimated values of Cdo and the Oswald efficiency, the aircraft drag polar was 
then compiled. 

Cd= c Do + Cl 2 /(it AR e) 

C d =0.018 + .0378 C L 2 

This plot (Figure 2.1-7) represents the total aircraft drag coefficient at the trimmed 
condition. The drag polar is of typical parabolic form, plotting the coefficient of drag 
versus the coefficient of lift for the entire aircraft. At the bottom of the 'drag bucket', 
where Cl=q, the parasite drag, Cqo, can be read directly from the graph. 


Drag Polar CD vs. CL 



Figure 2.1-7 


2 . 1-12 




2.2 


PROPULSION SYSTEM 


2.2.1 SYSTEM SELECTION 

In deciding upon the propulsion system best suited for this mission, we determined four 
important objectives which we needed our propulsion system to accomplish to satisfy the 
requirements of the mission. Therefore, the selection of the system was dictated by its 

ability to meet the following objectives: 

• provide necessary power for 3 full laps and return to landing area 

• maximize flight time 

• provide sufficient take-off power 

• be as lightweight as possible 

Initially, two methods of propulsion were considered. Because the system was required to 
be non-airbreathing and not emitting any mass (i.e. rocket, etc.), our options were limited 
to either rubber band power or an electric propulsion system. Analyzing the rubber band 
power alternative, we found that it would be difficult to obtain the necessary endurance to 
complete the entire course. Furthermore, the use of rubber bands would prevent us from 
having any speed control over the vehicle, since once the rubber bands are released, they 
cannot be stopped. However, in researching electric propulsion, it was found that an 
electric motor would be suitable to meet these requirements, as well as the objectives 
specified above. Thus, an electric propulsion system was selected for this mission. 

Next, the proper electric motor needed to be selected. After an initial analysis of the power 
required for the mission dictated that approximately 10 Watts of power are required during 
cruise (see section 3.3), it was decided that three motors would be analyzed: the ASTRO 
035, ASTRO 05 and ASTRO 15 motors. First, data was obtained for each of these 
motors. A propeller was selected with a known efficiency versus advance ratio relationship 
(namely the Zinger 10-6) to provide a starting point for the comparison of motors. Then, a 
computer program was written that provided power available data at various RPM settings 
(see Appendix A-l for listing and output of program). Next, the data for each motor was 
input into the program and the resulting power available data was analyzed. In looking at 
the ASTRO 035 motor system, it was found that this motor could provide up to 55 Watts 
of power once airborne and would have no problem keeping the aircraft in the air. 
However, in looking at the power required at take-off, it was found that the ASTRO 035 
would require approximately 96 feet to take-off, thus using almost 65% of the runway 
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simply for ground roll and exceeding the amount of runway that we wanted to use for take- 
off. Furthermore, this motor would require the use of a smaller propeller, thus diminishing 
the thrust available. Therefore, the ASTRO 035 was eliminated from consideration. 

The ASTRO 05 and ASTRO 15 were then analyzed in a similar fashion as the ASTRO 035. 
In looking at the two motors, it was found that each system would provide sufficient thrust 
both during cruise and in meeting the required take-off distance (see table 2.2-1). It was 
also found that the ASTRO 15 produced a significant amount more power than the ASTRO 
05. However, not only does the ASTRO 05 (6.5 oz.) weigh less than the ASTRO 15, but 
it also requires fewer batteries for operation and a much lighter total system weight. Since 
the additional power of the ASTRO 15 would be of no benefit because of the nature and 
goals of this mission, it was decided that the lighter-weight motor should be selected. 
Thus, the ASTRO 05 will be used to power the vehicle. 


Table 2.2-1 



MOTOR 

BATTERY 

SYSTEM 

Max. Power 

TAKE-OFF 

MOTOR 

WT. [oz.] 

WT. [oz.] 

WT. [oz.] 

Out [W] 

DIST. [ft] 

ASTRO 035 

4.5 

5.4 

11 

55 

96 

ASTRO 05 

6.5 

6.24 

13 

80 

31 

ASTRO 15 

7.5 

17.5 

25 

120 

17 


After deciding on the ASTRO 05, it was then necessary to select the type and number of 
batteries needed for the mission. Again, since one of our main objectives was to minimize 
the weight of the propulsion system, we wanted to use as few batteries as possible while 
supplying the motor with the necessary current. Due to the internal layout considerations 
of our RPV, it was also necessary to select batteries small in size. Furthermore, 
rechargeable batteries were chosen so as to keep the cost of batteries minimal if the aircraft 
is repeatedly flown. Using these guidelines, we utilized a TKSolver program (a 
programming tool found in the Notre Dame Aerospace Lab) which allowed us to vary the 
total battery voltage (number of batteries) and battery capacity while holding all other 
parameters constant (see Appendix A-2 for program ‘rules and output). Doing so, we 
found that a battery system of eight 500 MAH ('AA') Nickel Cadmium batteries would 
satisfy each of the objectives earlier defined. These cells will provide constant thrust 
throughout the flight, allow a motor current draw of 9.31 Amps and keep battery weight at 
only 6.24 oz. 
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2.2.2 PROPELLER SELECTION 


In order to achieve maximum efficiency at cruise velocity for the ASTRO 05 motor, a ten 
inch, two-bladed propeller with six inch pitch must be matched with the ASTRO motor. 
The final selection for the propeller resulted in the acquisition of the ZINGER 10-6. 

In trying to chose a propeller, the diameter is the most important parameter to consider 
since it effects every aspect of the performance of the aircraft. When conducting the 
analysis of several propellers, the diameter was constantly changed while other parameters 
such as battery voltage (9.6V), battery capacity (.5 amp-hr),weight (2.656 lbs), AR (12), 
and wing area (6 ft 2 ) are all held constant. The maximum diameter considered was twelve 
inches because this is the maximum allowable length to fit properly on the fuselage. The 
aircraft must have sufficient ground clearance for the propeller, and some reasonable 
landing gear lengths were chosen, which gave the maximum prop length. The minimum 
diameter considered was the eight inch because below this value the aircraft can not receive 
the necessary thrust to lift-off. 

A propeller-motor combination was needed that would achieve a short takeoff distance of 
approximately fifty feet with a large rate of climb, and have a flight time of over 120 
seconds. A quick takeoff is desired in order to reach the cruise height of fifteen feet 
quickly and allow ample time for the pilot to prepare for the first turn. A large rate of climb 
will allow the pilot to climb quickly after each turn in case any altitude is lost while 
executing the turns. 

The results of the Independent Trade Study, indicates that a ten inch propeller will deliver 
all of the needed abilities. The propeller-motor combination will takeoff in approximately 
31 feet, climb at 3.8 ft/s, and have a flight time of around 193 seconds (3.2 min.). 

Figure 2.2-1 shows the ZINGER 10-6 propeller efficiency vs advanced ratio curve. At 
cruise conditions the advanced ratio for the propeller is .277, which corresponds to a 
propeller efficiency of approximately .59. The curve was obtained from Dr. P. Dunn, who 
gave to his students during a design project in his class. 
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2.2.3 ENGINE CONTROL 

Once the aircraft has lifted off, the Drag-n-Fly will attempt to fly at a relatively constant 
altitude. During the straight portions of the course, we anticipate cruising at a gear RPM of 
7200. However, during turning flight, if we continue at the same RPM setting and flight 
speed, the aircraft will lose altitude. This occurs because as the aircraft turns, the vehicle 
will bank through the turns, changing the angle at which lift is acting and causing the 
vertical force on the wings to be diminished and consequendy, creating less lift Therefore, 
since we would like to maintain a constant altitude throughout the flight, an electronic speed 
control must be utilized. This speed control will allow the pilot to increase the RPM to 
approximately 8000, thus increasing the velocity of the vehicle and creating more lift over 
the wings. This additional lift will be used to maintain altitude through the turn. Once the 
turn is completed, the pilot will be able to decrease the RPM setting back to 7200 and slow 
the plane to its proper cruising velocity of 25 ft/s. 
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2.3 STABILITY AND CONTROL 

In the area of longitudinal stability and control the primary design parameter was the static 
margin A target value was set at 15% and the design entailed sizing the aircraft to ac e 
this goal. The other primary concern was trimming the aireraft a. the desired i lift _coeffic.e 
of 0 67 These areas will be presented in the nest two sections, 2.3.1 and 2.3.2. 
area of lateral stability and control, die primary concerns were with ^hievuigadequate ro 
control power with a combination of the rudder and the dihedral angle since the aircraf 
no ailerons. The equations governing dds condo, will be ^eloped and u^ o ^ 
vertical tail, dihedral angle, and rudder. With the resultmg large dthedral angle larg 
vertical mil. and a high wing, the roll and yaw stability of the aircraft was assumed more 
than adequate for the mission. Section 2.3.3 wiU present the analysis. 

2.3.1 STATIC MARGIN 

The primary goal in setting the longitudinal stability was achieving the desired static 
margin. The analysis of die static margin was performed in two separate phases. In the 
first phase die urge, sutic margin and center of gravity were set and the anidys.sdep* 
on sizing and locating the aircraf. components to attain the proper neutral pent, 
equation representing the neutral point of this aircraf, thus had to be derived and will be 
later in tins section. In equation form die definition used for die suuc margin was as 

follows. 

SM = Xnp/c - Xcg/c 

The second phase, or more detailed phase, used die same equation but die center of gravity 
and die neutral point were calculated and the resulting static margin determin 
using the equation above. This phase was used as die aircraf. was in die P^“ c " 0 " s “ 8 ' 
to check if die sutic margin was agreeing with the urge, value so changes could be made, 

possible. 

The static margin urge, value for this design was 15%. Typical aircraf, have a sutic 
margin between 5 and 10% bu, die urge, value was se, higher in this case to account for the 
different type of aircraft. The primary reason being the pilot is on the groun an nee 
extra subility for die RC aircraf, to allow a longer response time for die control. The urge 
of 15% was primarily used for sizing and placemen, of die tails. The urge, amertri «-» 
of gravity that was initially used was 35% of the mean aerodynamic chord thus urge, 
neutral point was 50% of the mean aerodynamic chord. 
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The detailed method used in determining the static margin was to estimate the neutral point 
and the center of gravity and take the difference between the two. The center of gravity 
was calculated from estimates of the sizes, weights, and locations of the internal layout and 
the components of the aircraft. This was discussed in the weights and structures portion of 
the proposal. The neutral point is defined as the furthest aft point where the center of 
gravity could be located without having a statically unstable aircraft. If the center of gravity 
was at the neutral point then the aircraft would be neutrally stable. The neutral point is 
found by setting the slope of the pitching moment equation equal to zero and solving for the 
center of gravity. The equation for the slope of the pitching moment curve is found by 
considering the contribution to the moment from each of the three major components of the 
aircraft; wing, tail, and fuselage. The following equation for Cma represents a typical 
aircraft 1 

Cm<x = Cm<xf - tiVhCloi (l-d£/dtt) + CLaw(Xcg/c _ X a c/c) 

For the RPV aircraft considered here this equation neglects one contribution. The 
contribution to Cma resulting from the moment created by the drag on the wing. This 
equation is valid for aircraft with the center of gravity a negligible vertical distance from the 
wing aerodynamic center. Since our wing was not at a negligible vertical distance from the 
center of gravity and the drag on the wing was high due to the low Reynolds number this 
contribution was considered. The equation for this contribution was approximated by 
considering moment as the drag times the vertical distance of the aerodynamic center minus 
the vertical location of the center of gravity. This equation then must be non- 
dimensionalized. 

Md = Drag ( Z ac -Zcg) 

Cmj) = Cdw ( Z ac /c - Zcg/c) 

From the definition of the coefficient drag 

Cdw = Cdo W + Cl ?/( n AR e) 
and the definition of the coefficient of lift 

Cl = Clow + Clow t*w 

At this point since we arc only interested in the contribution to the slope of the pitching 
moment curve, the terms independent of angle of attack, intercept terms, are neglected. 
Combining the previous equations and ignoring the portion that does not depend on angle 


1 Nelson, R.C., Flight Stability and Automatic Control, p.64 McGraw-Hill, St. Louis, 1989. 
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of attack produces the following approximation for the contribution to the coefficient of the 
pitching moment slope due to the drag on the wing. 

CmoD = ( 1/ (Jt AR e))( 2Clow Clow + Clow 2 iw )*( Zac/c - Zcg/c) 

A couple of approximations were made in calculations using this equation. First of all, the 
vertical distance of the aerodynamic center of the wing ( Zac) from the tip of the nose was 
an approximation of the location. Since the wing is a three panel polyhedral wing a large 
dihedral angle, Zac varies along this span. To make the equations easier to work with the 
value was averaged across the span. Also, in the above equation one term was dependant 
on the square of the wing angle of attack. For simplification of the calculations, this was 
approximated as the wing angle of attack times the incidence angle of the wing. Adding the 
previous equation to the previous equation for Cma produces the approximation of the 

pitching moment slope for our RPV. 

Cma = Cm<xf - qVHCLat (l~d£/da) + CLaw(Xcg/c-X a c/c) + Cm<xD 
Setting it equal to zero then solving for Xc g produces the following equation for the neutral 
point. 

Xnp = Xac/c + n Vh Qlcu/ ClccwO - de/da) - Cm<xf /C low - Cmao/CLaw 

The wing drag contribution comprises about 8% of the total coefficient of moment due to a 
change in angle of attack. The primary reason being the difference between the 
aerodynamic center and the center of gravity is about 6.5 inches or 7 5% of the chord. The 
above equation was used to size the aircraft components and orientations and also to 
produce the final approximations for the static margin as will be detailed in the next section. 


2.3.2 LONGITUDINAL STABILITY 

The primary goal of longitudinal stability was to achieve the targeted static margin, 15%. 
The other goal of the longitudinal stability was to trim at the necessary cruise lift coefficient 
of 0.67. The primary design parameters affected by these two conditions were the 
horizontal tail area, horizontal tail incidence angle, horizontal tail aspect ratio, and the 
distance to the tail aerodynamic center from the aircraft center of gravity. The distance to 
the tail had limitations though. The aircraft total length was set at 41 inches to minimize 
weight which limits the distance from the center of gravity to the tail to about 30 inches. 
The horizontal tail incidence angle was not effected by the static margin so it was primary 
set to trim the aircraft at the desired angle of attack. The wing incidence angle was not 
considered a stability design parameter as it was set at 5 degrees to produce the desired 
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cruise lift coefficient of 0.67. With the wing angle set, the fuselage will not have to fly at a 
positive angle of attack during cruise and extra fuselage drag can be avoided. 

The primary tool used in the analysis was the equation for the neutral point developed in the 
previous section. Since the target static margin was 15% and the target center of gravity 
location was 35% of the mean chord, the neutral point should be at 50%. Several of the 
parameters were varied and the design point that produced the static margin desired was 
used. 


2.3.2.1 HORIZONTAL TAIL SIZING, LOCATION, AND ORIENTATION 

Initially the horizontal tail area was set at 17% of the wing area of six square feet, or 1.02 
ft 2 , based on some RC aircraft rules of thumb. 2 Another area was used to note how the 
neutral point was affected by the different areas. This area was 20% of the wing area or 
1.2 square feet. Figure 2.3-1 demonstrates how the neutral point varies with the tail 
moment arm for the two different areas. 

Neutral Point as a Function of Horiz. Tail Area and Moment Arm 



Horizontal Tail 
Area (ft2) 

□ 1.02 

♦ - 1.2 


2 Lab Handout on rules of thumb for "Stability and Control" ,p.9, 1990. 
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^tailmomentannisthedistance 

gravity of the aircraft Note how the neutral potnt "roves back for ^ ^ ^ 

moment arm. The only point on this graph with a neu po 

moment arm of 2.0 fee. and a tail area of 1.02 square fee. as marked on die grap . 

for the two different moment arms. To yield P ^ ^ ^ 

aerodynamic chord of die wing, the 2.1 foot mome* J ^ necess ary is only 

required. When the tail area moment arm is mcreased 

0.83 ft 2 . 

Neutral Point Location as a Function of Horizontal Tail Are. 



Figure 2.3-2 


Tail Moment 
Arm (ft) 

a 2.5 

• 2.1 


. 0 ~ r. . the how the neutral point travels with respect to the 

One further plot, Figure 2.3.3 shows the how tne neu ^ ^ 

aspect ratio of the neutral point traveled 

... - ™^;~ “ 
aerodynamic chord Reynolds number. The flight Reynolds number 1.2 
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. „ .. . „ K _„, not v in5 Desoite this difference, no further 

tail local Reynolds number is about * 1U • Lres P uc 

problems associated with low Reynolds numbers were expected for this range. 

Neutral Point with respect to Tail Aspect Ratio 



Figure 2.3-3 

One further analysis that was done was to trim the aircraft at the desired lift coefficient of 
0 67 This was done by primarily altering the intercept term of the coefficient of the 
pitching moment curve. Again this included contributions ten the wing, tail, and fuselage 
but the primary controllable variable was the tail incidence angle. The results were 
produced by a simple computer program which used equations from the Stability and 
Control reference cited earlier. Oddly enough, the tail angle turned out to be about zero 
degrees (actually 0.6 degrees). This will make the construction of the tail easier as it wi 
not have to be tilted at any angle and can be just mounted relatively straight. The following 
figure 2.3.4 displays the moment coefficient as a function of the wing angle of attack or 
each of the components of the aircraft. These contributions were added together and the 
resulting Cm versus alpha relationship is noted as “aircraft” on the graph. The mm wing 
angle of attack of 6.5 degrees corresponds to a lift coefficient of 0.67 which is the desired 

trim lift coefficient. 
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Pitching Moment Coelf. tor Components of the Aircraft 



Figure 2.3-4 


The final values of the horizontal tail size, location, and incidence angle are summed up in 
the following Table. They were all based on the longitudinal stability goals of a static 
margin of 15% and trim CL of 0.67 although some values have been slighdy modified. 

Longitudinal Stability Summary 

Horizontal Tail Area 
Tail Moment Arm 
Tail Incidence Angle 
Horizontal Tail Aspect Ratio 

In the analysis it was assumed that the tail moment arm was constant. This was considered 
a fairly valid assumption for this type of aircraft since the center of gravity will move very 
little. The reason the center of gravity will move little is that no mass will be lost by the 
aircraft during the flight. No fuel will be burned and no mass dropped. 

The elevator effectiveness was also examined by considering its contribution to the pitching 
moment curve. The analysis demonstrated that a 15 degree elevator deflection would be 


1.00 ft 2 
2.3 ft 

0 degrees 

3.0 
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plenty to trim the aircraft from -2 to 12 degrees wing angle of attack. There is thus plenty 
of control power and maybe even too much. The aircraft will also be able to rotate to and 
trim at the 10 degree climb angle. Figure 2.3-5 displays the results. 

Pitching Moment vs. Lift Coefficient for Elevator Input 



Figure 2.3-5 


2.3.2.2 FINAL DESIGN 

In the final analysis the neutral point, using the equation given previously, turned out at 
55% of the mean aerodynamic chord. The center of gravity, estimated from simply 
balancing the plane, was at 45% of the mean aerodynamic chord. This meant a static 
margin of 10% as compared to the target of 15%. The reason the static margin was so 
small was because the center of gravity of the aircraft had slipped back during the 
construction of the aircraft. To alleviate the problem some of the internal layout of the 
aircraft were shifted forward and the wing shifted back to move the center of gravity 
forward. After some tests the center of gravity was moved forward even further by slightly 
weighting the front end of the aircraft. The final center of gravity location was at 30% of 
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the mean chord providing a static margin of 25® . This high static margin produced a very 
stable smooth flying aircraft 

2.3.3 ROLL AND DIRECTIONAL STABILITY AND CONTROL 

# 

The low Reynolds number RPV design that was proposed had to have a minimal weight in 
order to succeed. One step that was taken to reduce the weight was to discard the 
conventional ailerons. This eliminated the servos and control linkages for the two ailerons 
with the only penalty being a slightly larger tail and dihedral angle. The roll maneuvers 
then had to be performed with a rudder deflection. The main design parameters that affect 
this roll arc the dihedral angle and the rudder area. Other less significant parameters arc the 
lift coefficient, drag coefficient, and the vertical tail area. This section examines how the 
rudder induced roll can be approximated and how the design parameters effect this roll 
moment. In addition, the roll rate and how it was affected by the design parameters was 
studied. The process for the study was as follows. First, the governing equations were 
developed. Then, graphs were produced demonstrating the affect of the different variables. 
Next, the limits and constraints on the different variables were applied to produce a range 
of possible solutions. Finally, the most appropriate set of parameters was decided upon. 
In the analysis, it was assumed that the dihedral had the primary effect on the roll. All 
other affects such as the fuselage and the horizontal tail were considered negligible and 
were not quantified. In addition, as decided previously, the dihedral will be dispersed 
across the wing with a three-panel polyhedral. Thus this section includes the evaluation of 
the equivalent dihedral angle (EDA) and not its distribution across the wing. 

2.3.3.1 ANALYSIS TOOLS 

The theory behind the rudder controlled roll is that when the rudder is defected it creates a 
sideslip angle, p. As a result of the dihedral angle, the sideslip angle causes an effective 
change in the angle of attack of the wing. The leading half of the wing effectively sees an 
increase in the angle of attack and thus an increase in lift while the trailing wing sees a 
decreased angle of attack and decrease in lift The derivation few the expression quantifying 
this change in angle of attack as well as a diagram are given in Appendix D in the back of 

this proposal. The result is listed below. 

Aa = (3 sinr 

In this equation, T is the dihedral angle and Aa is the change in angle of attack. Using 
equation 1, an expression for the roll moment coefficient was derived. This derivation 
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began by assuming that the roll moment is equal to the change in lift on the wing, which is 
proportional to the change in angle of attack, times the moment arm. Page 1 of Appendix D 
includes this derivation that results in the equation given below. 

Cl = - 0.25 CL a sin(H P ( 1 ) 

The quantity CLa is the corrected lift curve slope of the wing. For a roll moment to be 
created a sideslip angle must be induced. A fixed rudder accomplishes this but does not 
produce a yaw moment as the RPV will trim at a non-zero sideslip angle. The equation for 
this relation is: 

Cn = 0 = Cnp P + Cn<i r dr (2) 

where dr is the rudder deflection. Equations for Cnp and Cndr were extracted from 
reference 1 and are listed below. 

Cnp = Cnp w f + ny Vy CL<xv (1 + dcr/dP) 

Cn^ = - ny Vy CL^xv T 

In equation 4, Cnp^f is an empirical factor representing the wing and fuselage contribution 

and was again estimated using a formula given in reference 1. Vv is defined as the vertical 
volume ratio ( lt*Sv/(S*b) ) and ny is the ratio of the vertical tail dynamic pressure to the 
ffeestream dynamic pressure. The term do/dp was approximated using an equation in 
reference 3 as cited below. 3 The parameter X is a factor which is related to the ratio of the 
rudder surface area to the vertical tail surface area. 

The only other tools used for the analysis were the LinAir program and the four part article 
“Dihedral” (reference 4). 4 LinAir was used to examine the affect of dihedral on the lift 
coefficient and the drag coefficient. The dihedral article was used to validate the data 
produced by the equations listed above and to produce a graph demonstrating the dihedral 
effect on roll rates. 

2.3.3.2 OBJECTIVES AND CONSTRAINTS 

The main objective of this study was to choose the dihedral angle and the rudder area. 
Merit was accredited to designs that produced an adequate roll moment which is defined in 
the next section. In deciding upon a dihedral angle, demerit was given to increased drag 


3 Nelson, R.C., Flight Stability and Automatic Control p. 71. McGraw-Hill. Inc., St. Louis. 1989. 

4 Beron-Rawdon, B. t “ Dihedral *. 4 part article, Model Aviation, Aug. ‘88 - Nov. ‘88. 
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Dihedral = 8 deg 
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each degree increase in dihedral is equivalent to about a 9% increase in the area ratio. The 
question at this point in the study was “Was this enough Cl to roll the aircraft?”. 

To determine an adequate range of Cl it was imagined that the aircraft had ailerons. The 
ailerons were sized using typical RC airplane data with reference to our aircraft. Then, 
using Equation 2.97 in reference 1 ( see footnote 1), Cl versus the aileron deflection was 
plotted for a maximum aileron deflection of +/-15 degrees. Figure 2.3-8 includes this plot. 
Based on this plot and the constraint for the aircraft that the rudder will not deflect more 
than 20 degrees, the necessary slope of the Cl versus rudder plot for our aircraft was 
determined to be 0.0024 per degree. Using this number the controlled variables could be 
varied so that the final design could be selected. 

The next set of graphs were the main tools in deciding on the final design for the dihedral 
and rudder. The slope of the Cl versus rudder input graph was set to produce adequate 
roll. Now all the parameters could be included on one graph to determine trends. The 
equation for Cl was of the form Cl=(A)dr where A was the slope. The equation for A is 
the following. 

A = 0.0024 1/deg = -1/4 CL aw sinT *(Cnd,/Cnp) 

All the parameters in this equation were defined previously. This equation could be solved 
for T which is called the rudder effectiveness parameter. The data for the plot was 
produced using the software Excel. The roll moment model data with the imaginary 
ailerons was typical data for a fairly acrobatic aircraft. This aircraft, requiring less roll 
control would not have to meet those standards. To correct this, it was assumed that only 
80% of the roll moment would be necessary (0.80*0.0024 1/deg). Figure 2.3-9 is the plot 
of the resulting data. Figure 2.3-9 indicates with the constraints that the dihedral angle 
would have to be something upwards of about 8 degrees and the rudder effectiveness 
parameter somewhere less than 0.81. As was noted previously, a one degree increase in 
the dihedral angle leads to a 9% decrease in the rudder effectiveness parameter. To find out 
which combination is the most feasible, more study had to be done. 

LinAir was used to determine if an increase in the dihedral angle to something above 8 
degrees would change the lift and drag significantly. Table 1 below gives the results for 
the modeled wing and horizontal tail at an angle of attack of 3 degrees. 
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Rudder Effectiveness vs. Dihedral Angle 
for Varying Vertical Tail Areas 


Vertical Tail 
Area 


-a 0.3 ft2 


0.5 ft2 
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Dihedral Angle (deg) 




Table 1 


Dihedral (deg) 

Cl 

C{) (wing and horiz. tail) 

8 

0.44521 

0.01887 

11 

0.43990 

0.01894 

14 

0.43295 

0.01883 


This data suggests that, for example, an increase in dihedral from 1 1 to 14 degrees would 
reduce the lift by 1.58% and increase the drag by 0.6%. This is not a significant amount if 
the dihedral amount needs to be increased. 

One other parameter that was affected by the dihedral amount was the roll rate. Based on 
some typical data given in footnote 4, the roll rate as a function of the dihedral angle was 
plotted for different yaw angles. Figure 2.3-10 is the plot. The roll rate that was desired 
for this aircraft was 20 degrees per second. This was based on a necessary bank angle of 
20 degrees, the cruise velocity of 25 ft/s, and a distance of 25 feet within which the bank 
angle is to be achieved. According to Figure 2.3-10, the only way the 20 degrees per 
second could be achieved for a yaw angle of 10 degrees is with a high dihedral angle 
(greater than 11 degrees). For a yaw angle of 12 degrees the desired roll rate could be 
achieved with a 9.5 degree dihedral. But a 12 degree yaw angle would need about a 24 
degree rudder deflection and the maximum deflection was set at 20 degrees previously. 
One other option would be to keep the yaw angle at 10 degrees and increase the speed of 
the aircraft as the roll is entered and then, decrease the speed and the rudder deflection as 
the roll is achieved ( only a 5 degree yaw angle is needed to hold the circular course). 

All calculations using the equations given in section 2.3.3. 1 were validated with graphs 
from the four part article “Dihedral”. The article supported the data produced and was in 
agreement to within a few percent. The LinAir drag and lift calculations were easily 
validated by calculating the change in the lift vector as it is tilted by the dihedral angle. 

Based on all the data and trends shown in the previous six graphs the following parameters 
were selected: a dihedral angle of 10.5 degrees: a rudder effectiveness of 0.68 which is 
equivalent to a rudder to tail area ratio of 50%: a vertical tail area of 0.50 ft 2 : a roll rate of 
18 degrees per second. The resulting rudder area will be 0.275 ft 2 . Note that these 
parameters give a roll performance slightly better than the 80% of the fictitious aileron 
performance (Figure 2.3-9). Because the roll rate is smaller than the desired value, the 
aircraft will have to increase speed entering the turns to achieve the proper roll rate. The 
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dihedral was not increased above 10.5 degrees because of the losses in lift and possible 
difficulties in construction and durability. The aircraft needs all the lift it can get to fly at 
slow velocities. The EDA of 10.5 degrees will be spread over a three-panel polyhedral. 
The outer section of the span will be 3.1 feet long on each end and have a local dihedral of 
12 degrees. The inner section will have an zero dihedral angle and be 1.2 feet on each side. 
With the large vertical tail and large rudder the directional stability and control was more 
than adequate for the mission. The large dihedral angle and high wing provided roll 
stability. A summary of the final design values for roll and directional stability and control 

are given in the Table below. 

Roll, Directional Stability and Control Summary 


Equivalent Dihedral Angle 
Vertical Tail Area 
Rudder Area 
Vertical Tail Aspect Ratio 


10.5 degrees 
0.50 ft 2 
0.25 ft 2 
2.0 


FINAL REMARKS 

The stability and control design of the aircraft is now complete. The design shows a very 
stable aircraft with plenty of control power. Based on how the technology demonstrator 
performed, the stability and control was adequate for the aircraft. The only necessary 
changes for future reproductions of this RPV would be to ensure the center of gravity is at 
30% of the mean aerodynamic chord. 
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STRUCTURAL ANALYSIS 
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Figure 2.4-1 


2.4.2 STRUCTURAL COMPONENTS 
2.4.2. 1 WING 
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The objective of the wing design is to produce a wing of minimum weight which would not 
fail under the maximum loading conditions during flight. This load condition is a load 
factor of 1.5 in a turn (48.2 0 in bank angle). Failure criterion is based on the axial stress 
due to bending at the root of the main spar, causing a shear failure. 

The design of the wing began with a Level 0 estimation. This estimation is based on prior 
database information for RPV’s. This information allowed the preliminary weight 
calculation and the preliminary design of the layout of the wing. It was from this analysis 
that it was found that the wing weight should not exceed 12 ounces and that the wing 
would consist of a single spar or multiple spars supporting the load. 

From the Level 0 estimation the Level 1 estimation takes a more in depth look at the wing 
design. This step involved the production of a short code on TKSolver! (a programming 
tool available in the Notre Dame Aerospace Laboratory). This code is used to calculate the 
bending moment at the root of the wing due to the lift on the wing. For this analysis, it 
was assumed that the spars of the wing carry all of the load on the wing as would a simple 
cantilevered beam. For this, half of the lift was placed at the midpoint of each side of the 
wing and the axial stress was found at the root based on the relation. 

Axial Stress = -Mx*v 
lx 

This equation is found in Beer and Johnston's Mechanics of Materials text p. 170. For 
greater detail on the analysis method and on the computer code used see Appendix B-l. 
The number of spars, rib spacing, t/c of the airfoil, leading edge and trailing edge size can 
be chosen and varied in this program. Some of these parameters do not enter the failure 
analysis of the program. For example the leading edge and trailing edge are assumed to not 
contribute to the strength of the wing, and hence are set to some reasonable size (based on 
database information, cross sectional area of leading edge = 0.0625 in 2 and trailing edge = 
0.25 in 2 ), and used only in the calculation of the total wing weight. The same is true of rib 
spacing which was set to 4” to maintain a smooth wing (spacing was set based on database 
information also), but has no effect on wing strength. Design of the wing depends on the 
factors that contribute to the spars ability to withstand bending. These factors are the 
number of spars, spar cap dimensions, and t/c of the airfoil. 

The purpose of the code is primarily to calculate the maximum axial stress in the spar to 
determine stress limits which must be met in order to avoid failure of the wing. Thus, a 
design figure of merit is the wing strength, as judged by a lower axial stress evident in the 
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Axial Stress vs. Height of Spar Cap: 
Lines of Constant Spar Cap Width 



Height of Spar Cap [inches] 


Analysis done with the following criterion set as constants: 

• Surface area of wing = 5 sq. ft 

• AR = 10 

• Span = 85.5 in 

• Chord = 8.5 in 

• t/c = 12% 

Figure 2.4-2 

The second part of the analysis involves relating spar cap cross sectional area to the total 
weight of the wing (See Figure 2.4-3). In this analysis wing surface area and t/c ratio are 
held constant. Since spar cap cross sectional area is obviously spar cap width multiplied 
times spar cap height, this graph can be used in conjunction with Figure 2.4-2. As one 
examines Figure 2.4-2 to select appropriate spar cap dimensions based on wing strength 
(resultant maximum axial stress), it is necessary to keep in mind the effect on weight that 
this will have, hence Figure 2.4-3 is used to see the relative effects on total wing weight. 
Total wing weight is based on calculations which include two spar caps, ribs spaced every 
four inches, a 0.25 in^ trailing edge, a 0.0625 in^ leading edge, spar webs 1/64" in 
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thickness, and a coating of Monokote. For these calculations spar caps were asssumed to 
be spruce, and everything else (other than the coating) was assumed to be balsa. 

Wing Weight vs. Spar Cap Area: 

Lines of Constrant Span 
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Spar Cap Area [sq. in] 


Analysis done with the following criterion set as constants: 

• Surface area of wing = 6 sq ft 

• t/c = 12% 


Figure 2.4-3 

At this point the third part of the analysis becomes useful. In Figure 2.4-4 spar cap 
dimensions and t/c ratio are held constant, while the effect of wing surface area on axial 
stress is examined. Lines of constant weight are added to maintain a reference forjudging 
any design. Wing design cannot be based solely upon strength and weight considerations. 
Performance considerations are the primary consideration in design of the wing. For 
propulsion reasons it would be desirable to have as large a wing area as possible in order to 
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Axial Stress [psi] 


Height of Spar Cap = 0.125 in 

Width of spar Cap =0.25 in 

Total Weight of Wing = 0.525 lb (8.4 oz.) 


Axial Stress vs. Surface Area of Wing 



Spans (b values) , and 
chords (c values) are in 
inches 

Weights (W values) are 
in pounds 

= Constraints 


Analysis done with the following criterion set as constants: 

• Width of spar cap = 0. 125 in 

• Height of spar cap = 0.25 in 

• t/c = 12% 


Figure 2.4-4 

After this Level 1 estimation was complete, the values determined for the wing were placed 
into Astros (a finite element analysis program) to verify their accuracy. This was used 
primarily as a check on what had been done previously, in order to ensure that the wing 
will not fail, and it was found that the wing was sufficiendy strong to avoid failure. 
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These pins are then analyzed as cantilevered beams. For more detail on the analysis of the 
empennage see Appendix B-3. 

2. 4. 2. 4 LANDING GEAR 

The landing gear of the vehicle is an important consideration. It is desired to have gear that 
will be strong enough to endure the landing loads but also have the ability to deflect and 
absorb some of the energy that will be transferred to the vehicle during landing. In sizing 
the landing gear for the RPV many parameters must be considered. The parameters that 
will be controlled in this particular analysis are modulus of elasticity, the length and the 
radius for the struts. The tire size was chosen based on research of the existing data and 
articles in model airplane magazines. The tire size was selected to be 2 in in diameter. 
Once these parameters have been solidified, the deflection of the gear due to the landing 
loads will be determined. 

In modeling the loads on the landing gear, the gear wire is assumed to be a circular 
cantilever beam with the load P acting at the end.(see Figure 2.4-6) The load is small 
enough that it may be assumed that the beam undergoes only elastic deformation due to the 
load. The give in the tire during the landing process has been assumed to be negligible 
compared to the deflection of the beam. The maximum deflection, x, due to this load 
occurs at the end of the beam and calculated by the following equation: x=P*L3/3*E*I, 
where E is the modulus of elasticity, L is the length of the beam and I is the moment of 
inertia of the beam cross section. This equation is found in Beer and Johnston's Mechanics 
of Materials text.The stiffness of the beam, K, must be determined and then used in the 
following equality to determine the deflection of the beam due to the landing load, 
l/2*m*v2=l/2*K*x2. In this equation, the mass of the RPV, m, and the landing velocity, 
v, are set by the structure and design requirements of the RPV and can not be varied. 
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A) Landing Loads 



B) Cantilever Model 
L 



Px 

Figure 2.4-6 


The stiffness of the beam varies with the radius, length and modulus of elasticity of the 
beam. A simple computer code was written to calculate the deflections due to loads varying 
from 10 to 100 pounds for varying radii of the beam with the length and modulus of 
elasticity as input. The code made use of the equations mentioned previously, (see 
Appendix B-4 for computer code and data) The data was then plotted and trends and 
restrictions were examined in order that a specific range of radii and lengths could be 
determined. 

The material of the beam was selected by holding the radius and length of the beam 
constant and calculating the stiffness from the load vs deflection curve. It is desired that a 
low stiffness will be used in order that the beam will bend and not remain stiff and transfer 
the load to the main wall in the fuselage. Observing Figure 2.4-7, it is obvious that the 
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Load vs Deflection 
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Figure 2.4-7 
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Figure 2.4-8 
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Load vs Deflection 



K=89.8 

K=46.0 

K=26.7 

K=16.7 


Figure 2.4-9 


Load vs Deflection 
Constant Length of 4 in 
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Load vs Deflection 
Constant Length of 6 in. 
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Load vs Deflection 
Conatant Length of 7 in. 
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2.4.2.S MATERIALS SELECTION 


In the design of each component in the RPV some type of analysis was performed in order 
to find the maximum stress during flight. From available data provided by the 
aerodynamics laboratory on some materials (See Table 2.4-1), the maximum stress can be 
compared with the maximum allowable stress of the materials, and an appropriate material 
can be chosen. From this, the RPV will be primarily made from balsa, except for some 
key structural pieces. These are the two spar caps, the wing hold down rods, and the 
vertical tail seat and pins, which will all be constructed from spruce. The three panels in 

the forward section of the plane are made of plywood due to the high loads in the payload 
area. 


= 9.0 

= 16.2 
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Table 2.4-1 


.. . tcptai . DENSITY [Ib/in3]_ 


Balsa wood 

Spruce 

Birch Plywood 


0.0058 

0.016 

0.023 


ALLOWABLE 
STRESS [psi] 


2.4.2.6 GENERAL CONSTRUCTION 
2.4.2.6.1 INTERNAL LAYOUT 

The internal components of die RPV consists of a motor, motor batteries, speed conuoller 
receiver, receiver batteries, two servos, and a payload <2” x 2” x 2 two ounce box). The 
design of the RPV retires tba, all of this weigh, be supported forward of the main ^ 
Jon of the fuselage. I, is also desired to tie these weights into iHewtngwhhw^ be 
supporting them during flight. For this design the motor ts placed as far forward 
possible, supported by a plywood fire wall on an end mount system. The -otor area t 
then connected to the main wing box area through a small box whtch ho ds t 
batteries. The main wing box is connected to the wing by rubber band hold down p . 
Two plywood panels a, each end of .he box add to its structural integrity. The w.ng box ,s 
also the area where most of the payload is contained. The 2 in cube is placed to the exneme 
from of the wing box, with the servos a. the exneme rear. The servos are s»PP«ted 
two bars across the width of the fuselage, while jus, behind the servos are the speed 
controller and receiver s, ached on top of one another, verrically. In between , he box 
payload and the speed connoUer lays the system baneries. velcroed down to the shelf. 

Figure 2.4-14) 

Top View of Fuselage - Internal Layout 


servos EJ controller 

_J l (system 
receiver 1 bat 


Figure 2.4-14 
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2. 4.2.6. 2 SCALE RPV DRAWING 

Since <he design team consists of amateur builders it was desired to keep the design of the 
“S as possible. As the design process developed and conttnued a « 

computer picture of the airplane was produced. I. was dnough the - 

picture that some construction difficulties were uienttfied and reso 
accomplishment of this task was to generate an accurate scale ptettue =of 
be used as a guide during actual construction (See Figures 2.4-5 a,b for plots of 
Hi ; exetc.se also allowed accurate placemen, of the shelf holdtng the sento 
"a, control rods will reach there respective control surfaces. also allows or 
determination of necessary connections between components, such as the vetttca. tatl sea . 
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2.5 WEIGHT AND CENTER OF GRAVITY ESTIMATION 


Much of the initial work done by the weights and structures division was to estimate the 
weight and center of gravity for our RPV. Due to the fact that the weight of the RPV was a 
primary concern throughout the design process, it was necessary to spend the effort and 
time in keeping the RPV systems and structure to a minimal weight to increase the 
possibility of a successful mission. The initial step was to investigate the existing data 
base. Once this was done and a "preliminary" estimation of the weight of the RPV was 
completed, the center of gravity could be calculated. It was desired that this center of 
gravity fall within the range of 25% to 40% of the chord of the wing in order to obtain the 
desired stability and control for the RPV. Of course as the design progressed and different 
aspects of the vehicle changed, these values also changed. The design group set initial 
goals for each particular division, for example structures and weights, propulsion, and 
stability and control, to meet. Once the data base was reviewed, the design weight for the 
RPV was set at 3 lbs (48 oz). The structures and weights group felt that if this design goal 
was met, as well as the other design goals, the mission would be successful. 

2.5.1 COMPONENT WEIGHTS 

As stated previously, the first step that was necessary in this particular analysis was to 
examine the existing data base. Speaking to our 'resident expert' on RPV systems, Mr. 

Joe Mergen, a good estimation for the sizes, material and weights for RPV components 
was obtained. From this review of the data and making estimates based on this data and 
the desired size of our RPV, a "preliminary" estimation of the weight was made; level one 
of the weight estimation. 

Level two of the weight estimation for the vehicle consisted of obtaining more concrete 
values for the weights of the components of the RPV, such as the motor, servos, radio 
receiver and batteries. Once these have been determined, the size and weight for the 
structure of the vehicle was estimated. Again the existing data base was consulted, but 
values were also obtained from the manufacturers of these components once these 
components were selected. The propulsion team selected a certain motor and the required 
batteries to operate this motor. The number of servos was determined by the stability and 
control division. The initial estimation for the weight of the propulsion system was 25.0 
oz, which includes the battery pack and motor. The remaining avionics were initially 
estimated to weigh 18.5 oz. 
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Once the weights for the propulsion system, control system and remaining avionics were 
estimated, the weight for the RPV structure was estimated. From the des.gn go s 
determined early in die design ptocess, it was noted that the total wetght for the RPV was 
Leal factor in completing the mission successfully. The weigh, and sm.ctures drvrs.on 
of the group decided that in order to keep this weight a minimum, the av.on.es and 
propulsion system would be selected and a minimal structure would be constructed around 
it speaking wth Mr. Merged and examining RPVs burl, in previous years, the 
materials were tentatively selected and the general structure of the RPV detemuned. 
initial estimation for the wing weight was 12.0 or and the weight for the fuselage / 
empennage section was 12.0 oz. This brought our initial estimation to a value of .3 

(68.8 oz) for the total weight of the RPV. 

The wing section was then studied in greater detail and was designed to be a wing of 
minimal weight that could withstand normal flights loads. This- design process .*» 
discussed previously and resulted in a second esdmadon of the wmg we.ght of 8 0 oz. As 
d,e design of the vehicle proceeded this value changed. The final esdmanon for the wetght 
ofd.ewingis8.4oz. The fuselage / empennage section was then exammed .n greater 
detail It was decided by the weights and structures division of the group that the fuse g 
would be of a truss configuration in order to keep the weight down by reducing the nee 
for solid sections in the fuselage structure. Solid panels would be tncluded i .n the mam 
wing box section, where the propulsion system and avionics wtll be located for strong 
reasons and the empennage was designed to be a minimal strocture that could withstand the 
load incurred during the mission. The sizes and marerials for the members of the mtss and 
other portions of the fuselage were determined by examining model amplane magazmes and 
consulting Mr. Mergen. Once this stracture of minimal weight was destgned .t was 
analyzed to insure that failure would not occur during flight and a value for e we.g t o 
this fuselage / empennage rection was estimated to be 7.2 oz. Again, as the, design 
RPV proceeded the weight of tile fuselage / empennage section changed. The final 
estimation for tile weigh, of this section is 8.4 oz. The landing gear was also ^d.ed .n 
greater detail and resulted in an estimation of a weigh, of 3.5 oz. As was stated prevtously, 
me structure would be changing as me design progressed and thus the wcghts of the* 
components would be chanpng. The following table. Table 2.5-1, shows me final values 
for the weight and weight percents for the components of the RPV. 
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TABLE 2.5-1 


COMPONENT 

WEIGHT (oz) 

WEIGHT PERCENT 

Avionics 

8.0 

18.3 

Payload 

2.0 

4.5 

Motor 

12.7 

29.1 

Propeller 

0.7 

1.6 

Wing 

8.4 

19.2 

Fuselage/Emp 

8.4 

19.2 

Landing Gear 

3.5 

8.1 

TOTAL WEIGHT 

43.7 oz 

DESIGN WEIGHT 

48.0 oz 


2.5.2 CENTER OF GRAVITY ESTIMATION 

Once the component weights were estimated, the center of gravity could be located. This 
calculation was done several times throughout the design process, due to the fact that it was 
a necessary piece of information for the stability and control division. The weight and 
location for the individual components must be known in order to calculate the center of 
gravity. Therefore it is essential to study the internal layout for the systems and avionics in 
order to position them correctly. The systems were positioned and mounted in such a 
manner that they would be easily accessible and push the center of gravity as far forward as 
possible. This was desired due to the short nose of the vehicle. Due to this short nose, it 
was a concern of the weights and structures division that the center of gravity would be too 
far rearward. 

After the initial estimates of weights were refined through the design process, the center of 
gravity was found using locations and weight figures for each component of the RPV. 
These component weight figures are more accurate than the initial weight figures, since the 
initial figures were based largely on preliminary estimates. A table summarizing the 
majority of the RPV's component weights and locations is given below, Table 2.5-2. The 
component center of gravity location is referenced from the nose of the RPV. (See Figure 
2.5-1) 
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TABLE 2.5-2 


COMPONENTS 


CENTER OF GRAVITY 
i nr ATION (in) . 



Fuselage (main truss) __ 


Landing Gear 


Horizontal Tail 


Vertical Tail 


Wing Box 


KT/^cf. r mot or casing) 

- T i 37 1 J,M — 

Nose (batteries) I - 

Xcg from nose = 9-4 in - 
RPV it is expected to fall near the X axis. 
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Side View Fuselage - Component Location 


speed 



Figure 2.5-1 
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3 PERFORMANCE ESTIMATION 

3.1 TAKE-OFF AND LANDING ESTIMATES 


Take-off distance for the aircraft is constrained to only 150 feet. This is due to the fact that 
the flight will occur inside the LOM US Sports Complex (Notre Dame, IN) and space is 
limited. With the proposed propeller-motor combination, a take-off distance of 31 feet is 
achieved. This value is based on equations found in Introduction To Flight (p. 309-311) 
by Anderson. 

take-off distance 1.44 * W 2 

g * P * s * q, ma **T[D+ p r(W-L)] 

When using the equations, ground effect was neglected because the effects are beneficial 
and a worse case scenario was used to insure that the aircraft would lift-off. Also, a 
ground friction coefficient of .04 was assumed for the astro turf found in the LOM US 
center. 

Once the aircraft completes the required three laps is must prepare itself to approach the 
runway and land. Upon the wheels touching the ground, the aircraft will rely strictly upon 
ground friction as a means to slow the aircraft to a halt. Using the equations found in 
Anderson’s Introduction To Right (p. 313-315), the total ground landing distance is 133 
feet. 


landing distance 1.69 * W 2 


g*p*S*C Umax *(D*p r *W) 


This distance is the actual distance it will take the aircraft to stop once the wheels touch the 
ground. A ground friction coefficient of .04 was once again assumed for the astro turf. 


3-1 



3.2 RANGE AND ENDURANCE 


In evaluating the range and endurance of our RPV, it must be noted that these parameters 
are of importance because the vehicle must be able to complete the full three laps as slowly 
as possible and return to the landing area. However, excess range and endurance 
capabilities are of no extra benefit and, in fact, could be detrimental to the goals of the 
mission if these “benefits” caused an increase in the RPV’s weight Thus, it was necessary 
to try maximizing endurance and range by using the aerodynamic relationships and 
conditions in R. C. Nelson s Right Mechanics text (pp. 4.16-4.17). However, since it is 
impossible to simultaneously maximize range and endurance, it was decided that endurance 
would be maximized since this was one of the design objectives specified in section 1.1. 
This was accomplished by trying to minimize power (thus maximizing Cl^/Cq). 


To actually estimate the range and endurance, a TKSolver program was again utilized to 
determine the sensitivities inherent to these performance characteristics. As flight speed 
and propeller diameter were varied, it was found that the endurance proved to be very 
sensitive to a decrease in propeller diameter. This is evidenced by the fact that a three inch 
decrease in diameter increases endurance by nearly 1 10 seconds (see figure 3.2-1). 

Endurance vs. Prop. Diameter For Varying Flight Velocities 



This corresponds to the fact that as the propeller diameter decreases, the motor does not 
have to work as hard to drive the propeller. Thus, motor current draw decreases and 
endurance increases. As a result, it is estimated that the Drag-n-Ry can stay in the air for 
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„ t0 ,93 seconds on a single charge (Ml battery capacty) at an airspeed of 25 ft/s and a 
propeller diameter of 10 inches. 

A similar analysis was done to determine the ^^^^^^^g^^anwter'lseefi^ue 
,ha. the general trend is for range to tncrease with decree g P determined 

tot the RPV must travel approximately 2800 feet from distance in 

1, stops after landing. Consequently, our atrcraf. must ^ fl y l p to 4831 fee. a. 
order to complete to mission. Thus, it is estimated “ ^ gh the ^ is 

an airspeed of 25 ft/s and with a propeller excess -ge did no, 

substantially greater than tot ' 1 earher) and provides a comfortable 

adversely effect the weight of the RFV tas 

margin of safety. 



3.3 


POWER REQUIRED AND POWER AVAILABLE 


In analyzing to power required for to RPV, a number 

Parameters such as aspect ratio, annual power required 

parasite drag all had to be estimate as written which accepted 

curve. Once these parameters were set. a computer program was w^ „ 

these values as input and ton varied airspeed over a spectfied range (see App 
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listing and output). Using this program it can be seen from figure 3.3-1 that at our cruise 
velocity of 25 ft/s, 7.83 Watts of power are required to power the RPV. Furthermore, at 
25 ft/s we are a safe distance from the minimum power required, thus eliminating any 
concerns about “bottoming out” the power required curve, but also maximizing endurance 
as much as possible. Bottoming out the power required curve is a problem because at 
minimum power required, the vehicle’s speed is not very stable and has a tendency to vary. 
Therefore, for effective speed stability, we did not want to fly at precisely the minimum 
power required. At take-off, our calculations determined that 46.2 Watts of power are 
required for a worse case situation (no ground effects contribution to lift). Thus, it was 
necessary to ensure that we had at least 46.2 Watts at all times. Figure 3.3-2 demonstrates 

that at our settmg of 7200 RPM, our power available is safely more than our power 
required. 


POWER REQUIRED CURVE 
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POWER REQUIRED AND POWER AVAILABLE AT 720* RPM 



3.4 


CLIMB AND GLIDE PERFORMANCE 


The aircraft must complete and extra 

for landing at the end of the mission. g tance covere d in the air to finally reach 

sixteen degrees and **^*^££ descent is begun immediately after the final 
the ground is 52 feet ( dist - h / tan Y)- 11 111 com e tQ a 

turn, the aircraft will have enough distance to ^ feet ofrunway 

complete stop. This is possible because o te ^ stop ^ aircraft once on the 

distance is ample distance for t c air ^ ra _ ad ds t0 enough room for 

ground is 133 feet and the distance in the air is 52 feet, whicn 

the aircraft to glide and stop completely. 

™ aircraft's sea .eve, ctaib and 

shown in figure 3.4-1. The rate of climb at take ott PP 
maximum rate of climb is 4.1 ft/s a, a flight velocity of 23 ft/s. 
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RATE OF CLIMB AND DESCENT 
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4 technology demonstrator 


4.1 FINAL CONFIGURATIO ^ to the original plans. It must be 

The final flight aircraft rolled out with a n ^*ro ^Lnt to building accurately to the 
noted that computer aided drawings were P necessity. Holes were 

oHsin, design. CHan.es of an ounce for each o, 

cut in the rear spar of the mam wing torsional strength. 

Table 4.1-1 

^ i t Measured Demonstrator weights 


COMPONENT 

Fuselage . , , 

(Avionics-Motor included) 

Full Wing 

Horizontal Tail 

Vertical Tail 

Landing Gear 


WEIGHT (grams) 
800.0 

210.0 

49.5 

25.2 

57.0 


WEIGHT (oz) 

28.22 

7.41 

1.75 

0.89 

2.0 


Total (by component) 
Total (final) 


40.27 oz (2.52 lb) 
48.8 oz (3.05 lb) 


, tll _ final weight. There are several 
The sum of ft, component weights does no. was appUed . Lead 

reasons for this. The components were we.g Also, a brass 

weight weighing 4.0 oz was added to bafcnce « * “ * J 
spinner (weigh, unknown) was put on to place to pt»pell« 
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4.2 TECHNICAL DEMONSTRATOR FLIGHT TEST 

With fabrication con-pice, the High. Test program had as "m 

Z, K s< had specific objectives to accomptish and a flight pian to meet dtese objecuves. 

Due to the unpredictability of the aircraft performance, oniy the design team members and 
rr^ af-ed on the fteid. Ali spectators watched horn behind the track, or from 

the bleachers. 

Tuesday, May 1: Taxi tests 


Objectives: 

1. Validate ability of aircraft to perform takeoff roll. 

2. Determine the basic handling characteristics of the aircraft 


on the ground. 


Flight Plan: 

1. Run aircraft up Field. 


2. Quick take-off and short glide if possible. 


Attempt 1 

Center of Gravity placement. 

. Check shows C.G. at approx. 40 % chord. 

• C.G. Placement should be at main spar. 

• Conclusion: Placed too far rearward. 

. Solution: ballast the nose with coins and keys and move engme forward 
approximately one inch. 


System checkout: 

• Propulsion okay. 

• Controls range check okay. 

• Landing gear check: 

• gear bowed outward. 

• questionable stiffness. 

• provisional okay. 
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Performance: 

• Attempt to taxi unsuccessful. 

• Propulsion could not overcome 


roughness of Astroturf due 


to bowing of landing 


gear. 

• Solution: 


Attach two more pieces of wire to 


the fuselage to make landing gear 


more rigid. 


Attempt^ 

. C.G. placement : soil too far forward, though improved. 
. Systems check okay. 

• Landing gear okay. 

Performance: 

• TO distance of 20 yards. 

• Climb to 8-12 feet. 

• Power cut. 

• Glide to landing. 


Conclusions: 

1. Objectives met. 

2. Demonstrated capability to taxi, takeoff, and glide. 
Wednesday, May 2 : modifications 

. Lead ballast bringing C.G. to correct position painted and taped to nose. 

• Wing moved as far back as possible. 

• Additional wire strung to landing gear. 

Thursday, May 3 : test 

Attempt 


Objective: 

Repeat successful test of Tuesday. 
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Systems Checkout: 

• C.G. and all systems check out. Some concern over slight wing lean. 
Performance: 

• Takeoff Distance of 20 yards. 

• Reached altitude of 10 ft. 

• Engine cut after flying in straight line. 

Conclusions: 

1. Demonstrated: 

• Repeated taxi capability. 

• Repeated takeoff capability. 

• Extended glide 10 yards. 

• Flight in straight line. 

• Dynamic stability. 

2. An unqualified success! 

Friday, May 4: Final Flight Test 

Objectives: 

1. Validate turning capability (right and left turns). 

2. Perform 3 figure eights as required. 

Attempt 1 

• No further changes to aircraft from Thursday. 

• Successfully completed the proposed flight path. A total of six figure eight laps were 
made. However, this flight was run at a velocity much faster than desired. 

• All systems checked out okay. 

Performance: 

• Total flight time: 0:01:53 hours. 

• Pilot started plane slowly for approx 10 yards. 

• Takeoff distance was approximately 20-25 yards (total 35 yards). 

• Flew six full figure eights. 

• Maintained level, controlled flight; no imbalance in wings. 
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• Full turn approximately 60-7 5% of field width (40 yards). 

. After the final figure eight lap, the plane penned a mm to return to level fhgh, 

and landed before batteries emptied. 


Conclusions: 

1. Both goals accomplished. 

2. Full flight requirements accomplished. 

3. Established: 

• Turning capability. 

• Dynamic stability. 

. Endurance to complete mission without power failure. 

Attempt! 

Objectives: 

1. Quantify flight speed and altitude. 

2. Fly again to establish repeatability. 


Plan: 

1. Takeoff. 

2. Fly oval around field. 

3. Fly figure eight. 


Performance : 

• Total flight time 0:01:35 hours. 

• Cruise velocity : 24 ft/s. 


. Re = 99,000. 

• Take-off distance of 75 ft 

. Hew to level with nets approx 30-35 ft.while flying in an oval path. 
. I .Irin.de to 20 feet to prepare for landing approach. 


Conclusions: 

1. Both goals accomplished. 

2. Established repeat performance. 
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5 feasibility of high altitude 

MISSION 

The basic objective of this low altitude design study was to determine some ot the problems 
Z involved in low Reynolds number flight. Iw,w Re fltght has tmp catto. 

for the design of aircraft and estimation of the loads and drag forces u will encounter Ove 
the course of the s.udv, many of these implications were discovered; they came out w te 
normal techntques of esttmation of parameters were applied to situations that are out o 
scope of the ordtnary. These impl.ca.ions are also valid for low Reynolds fltght a, t 

altitudes. 

A high altitude tntssion could encounter low Reynolds flight, both because of lower speed 
IXs in air denstty . For example, the dens,, of air a. 10,«0 leetis t tu of s 
level Assuming the viscosity is constant, this would drop the Reynolds number by . 
A, 75 000 feet, the density is 5% of sea level. The station keeptng nature of the mtsston 
would require low speeds, as tt ts difficult to fly in a tight area at htgh speed 

One of the problems of low Reynolds flight comes at the level of data sources. 

Commotdy, the espenment* data cotrnng ou, of wind tunnel research oes not „tend tnto 
the range of Re=100,000 unless that level or lower is the target o esu 
much NACA data are taken at an Re ten times greater than this, Re=l * 10 an 

aerodynamics, the implications come with estimadon of the atrcraf, Itf, charac, ensue. 
Infinite wing airfoils can deliver c<*ffic,en,s of lift in etccess of 1 .0. However, » - "y, 
the true wing lift coefficient can be much smaller, in the ranger of 0.5 - 0.8. A . rc rat, 
designed for higher speeds (htgher Reynolds number) can overcome thts loss by Ilyu g 
a faster velocity. At low Re, the lifting area of the aircraft must be larger to otfse. the 
corresponding loss of lift associated with low speeds or low denstty. 

Drag ts another area affected by low Re flight. The data base is one of the problems here. 
The key problem can be seen from lookrng at the drag polar. As the aircrat, goes to lower 
coefficients of lift, to the left of bottom of the drag bucket, the drag actually goes up. 
Increased drag can then be a problem. 

Propulsion is affected also by lower velocities and densides. At low Re, the required 


power, a direct function of drag, also goes up. Power required is also inversely 
proportional to the square root of the density. As density decreases, the power rises. 
Therefore, low Reynolds flight increases the demands on the propulsion system. 

Stability and control is particularly sensitive to low Re flight. Control surfaces must be 
enlarged to account for the loss in effectiveness due to the decrease in air flowing by the 
surfaces. The surfaces cannot effect as much influence as at higher speeds or lower 
densities. 

The design study has been useful for exploring the difficulties of design and analysis of a 
low Reynolds flight vehicle. The successful flight of the Technology Demonstrator shows 
that these difficulties can be overcome. 



Appendix A-l 

Power Available Analysis and Computer Code 


To determine the power available for each particular motor, it first was necessary o vary 
velocity and RPM using a DO loop in the code. After a velocity was chosen, an RPM 
setting was then selected. Knowing these parameters, the advance ratio could 
determined by assuming a propeller diameter: 


Advance ratio = V/nD 


where D is 10 inches because for this progrom. the Zinger 10-6 was arbitrarily selected. 
Next since the velocity vs. advance ratio data was available for the Zinger 10-6,.he 
efficiency was detetmined for the propeller at this velocity and RPM setting. The gear 
power was then determined by using die dam provided by ASTRO for each motor, creating 
a plot of engine RPM vs. gear power and then determining the equation of that curve 
this manner, it was then possible to approximate gear power solely as a funcuon of RP . 
Finally, die actual power available for this velocity and RPM was found by multiplying the 
gear power and propeller efficiency. The program then selected new values of velocity an 

RPM and the process was repeated. 



SALFORD UNIVERSITY 
COMPILER OPTIONS: 


<USER2>S424049795>STUDENT>PR F77 


FTN7 7 VER. 22QS _ 

<jad momad rHECK MOB I G LQGL DYNM NOOFFSEl L»' 
EjB^E^SJ^SlSE NOIMPURE 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 
CO 19 
0020 
0021 . 01 
0022 . 01 

0023. 01 

0024. 01 

0025. 01 

0026. 01 
0027. 01 
0028. 01 

0029. 01 

0030. 01 

0031. 01 

0032. 01 

0033. 01 

0034. 01 

0035. 01 

0036. 01 

0037. 01 

0038. 01 

0039 

0040 

0041 

0042 

0043. 01 

0044. 02 

0045. 02 

0046. 02 

0047. 02 

0048 02 

0049 02 

0050. 02 

0051. 02 

0052. 02 

0053. 01 

0054 

0055 

0056 

WARNING - 


C 

C 

C 

c 

c 

c 


IWR 


REAL X { 31 > , Y(31. 1 > 

WRITE (1.*) 'ENTER 
READ(1.*> IWR __ A - 0£rrT 

WRITE <1,*> '| N I E g ?nO ECT 
read(i,*> ar.s. e.w.cdo 
WRITE (1,*) 'ENTER DIAMETER 
READ ( 1. * > DI 


RATIO. S (FT**2 ) . E. W (LBS), CDO 
OF PROP IN FT' 


WRITE! IWR, *> 
WRITE! IWR, *) 
WRITE! IWR, *> 
WRITE! IWR, *> 
WRITE! IWR, *) 
WRITE! IWR, *> 


'ASPECT RATIO*', AR 
'PLANFORM AREA (ft**2>* 
'EFFICIENCY* ', E 
'WEIGHT ! lbs )* « W 
'CDO* CDO 


V! f t/s ) 
D! lb ) ' 


PR!W) 


10 


WR I TE ! I WR , * ) 

WRITE ( IWR, *> ' 

CD 

DO 1 10 V=10, 40 

X ! I ) =V ,,, _ . 

CL=W/!. 00119* V**2. )*S) 

CD=CDO+ ( CL**2 ) / ( 3. 14159*E AR > 

D=0 00119*CD*V**2. *S 
PRI=D*V 

Y! I, 1 )=PR 1*1 • 356 
PR*PR 1*1 356 

WR I TE^T IWR < * ■ V.PR.CL, CD, D 
CONTINUE 

I OPT*— Oil 
ND-31 

paLL TPLOT ! IOPT, X, Y, 31 , ND, NF ) ... . . , 'POWER 

gJtt T*ITLeV 'PRELIM I NARY^POWER REQUIRED' CURVE ' ) 


CL 


V! ft/s > 


500 


WRITE (IWR,*) 

WRITE ( IWR, *) ' 

EAL !W) PG 

DO 98 V=10, 40 
DO 99 N=4700, 9000, 

AN=N/60 

EP*+ / 07807+2. 5764*AJ-2 7452*AJ**2 +. 
r PG=-80 062+ . 11418*N-1. 

C PG*-88 153+. C8739*N-8. 5894E-6*N**2. 

PA=PG*EP 

PAREAL=PA*1. 356 oao _ a| pr 
WRITE! IWR, * > V, N, EP, PAREAL, PG 
99 CONTINUE 
98 CONTINUE 
77 CONTINUE 
STOP 

_abel 77 E has not been referenced 


RPM 


REQUIRED (watts 


EP 


86713*AJ**3. 


END OF COMPILATION CLOCKED 


457 SECONDS 
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Appendix A-2 

Battery Voltage/Current Analysis and TKSolver 

In deciding «*>n the desired baiter, voltage and ' 

code was implemented which was provided by ^ ^ other 

This code allowed for quannnes such asve 8 ^ ba(tery capacity cou id 

pertinent parameters to be set. ce - effects of these variations 

be varied while keeping all other parameters «-«£££ „ umbCT rf battcri es) 
could be studied. While varying .he battery ™ Uag ' ^use a motor 
audbanerycapacity ^£^^2^***** 

current draw in excess of P a , s0 had to be monitored so as to ensure 

speed control. Furthermore, range ft was ^ parameters (range 

the vehicle could meet the o jective number and type of batteries 

and endurance) that were the driving factors m selecnng die number and typ 

for the Drag-n-fly. 


S Rule 

* 0“ .5*rho*vel A 2 

* Cd-Cdo*CI A 2/(PI 0*eff*RR) 

* CI-(n*U)/(Q*S) 

' req*Q*S*Cd*ve I 

* ROC*(Pauoi l-Preq)/U 

* u"uset-Kb*i 

* »otrp» - (y-i*Ra)/Ky 

* proprps*«otrp«/(60*gr) 

* >vel/(proppp3*propd) 

* CT-Ct(J) 

* CP-Cp(J) 

* eto-Ct(J)*J/Cp(J) 

* Pavai l ,, eta*Cp(J)*pho*ppoppp3 A 3*ppopd A 5 

* Cp( J)*pho*ppoppp« A 3*ppopd A 5-( (Kt/Kv)*(v* i - 1* !*Ra)* . 0005454- f I o«(«o 

* f ltl«o*batcop/l 

* pange - ye I * f 1 1 1 «e*3600 



Input 

Note 

Output Unit 

Coaaent 


Q 

.713125 pef 

dynaaic pressur 

.002378 

rho 

9 1 ug/ f 1 3 

air density 

25 

vel 

ft/sec 

air speed 


Cd 

.05176516 

a/c drag coefficient 

.0-1 

Cdo 


zero lift drag coefficient 


Cl 

.59568265 

a/c 1 i ft coefficient 

.8 

eff 


efficiency factor 

12 

RR 


aspect ratio 

1 

n 


load factor 

2.656 

U 

lb 

a/c aeight 

6 

S 

ft-ft 

aing area 


Preq 

7.822981 U 

a/c poeer required - level f 


ROC 

3.7636216 ft/s 

rate of cl iab 


Pava i 1 

21.375295 U 

poaer available froa propel 1 


v 

8.6116061 volt 

araature voltage 

9.6 

vset 

volt 

battery voltage 

.1058 

Kb 


battery constant 


1 

9.3137391 aap 

aotor current draa 


»otrp» 

16111 .561 rpa 

aotor speed (rpa) 

.06 

Ra 

oh« 

araature resistance 

.0005 

Kv 

vo It /rpa 

aotor speed constant 


proprps 

7290.3005 rpa 

propeller speed (rps) 

2.21 

9^ 


gear ratio 


J 

.2170022 

propeller advance ratio 

.833 

propd 

ft 

propel ler diaeeter 


eta 

.58833081 

prope 1 1 er ef f I c i ency 

.62 

Kt 

in-oz/aap aotor torque constant 

.95 

greff 


gear efficiency 


fit lie 

193.26288 sec 

flight t i ae 

.5 

bat cap 

aap-hr 

battery capacity 


range 

1831.5719 ft 

range 


CT 

.03730893 



CP 

.01566362 




Appendix A-3 

Power Required Analysis and Code 

To determine the power required curve for the Drag-n-Fly, a computer code was written 
that allowed for the Aspect ratio, planform area, Oswald Efficiency factor, weight and C^ Q 

to be input into the program by the user. A DO loop was written to vary velocity from 10 
ft/s to 40 ft/s and for each run through this loop, the lift coefficient was determined by 
assuming steady level flight (L=W): 


C L = W/(Q*S) 

The drag coefficient was then determined by adding Q )0 and the induced drag: 

C D = Q,o + C L 2 / JteAR 


Next, the total drag was calculated from the drag coefficient above: 


D = C D *Q*S 


Finally, the power required is simply the product of the drag and velocity: 


Preq=D*V 


Using this process for each run through the code, the power required curve seen in Section 
3.3 was produced. 



SALFORD UNIVERSITY" FTN77 VER. 228S <USER2>S424049795>STUDENT>PR F77 


COMPILER OPTIONS: 


15: 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
001 1 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 . 
0022 . 

0023. 

0024. 

0025. 

0026. 

0027. 

0028. 
0029. 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 
WARNING 

0041 

0042 

0043 

0044. 

0045. 

0046. 

0047. 

0048. 

0049. 

0050. 

0051. 

0052. 

0053. 

0054. 

0055 

0056 

0057 

0058 
WARNING 


N0FRN N gpJ N & L 8&REC SoSILENT E NO_OPTimSE G NOIMPURE OOFFSET L ®° NC 


REAL X ( 31 ) , Y(31, 1 ) 

oSIIf, (1 '*> 'ENTER IWR ' 
READ< 1, *) IWR 

WgJJf 'ENTER ASPECT 

READC1, *> AR, S, E, W, CDO 
WRI If, (1 ' *> 'ENTER DIAMETER 
READ ( 1 , *) DI 


RATIO, S (FT**2),E, W (LBS), CDO 


OF PROP IN FT' 


WRlTEdWR, *) 
WR I TEC IWR, *) 
WR I TE ( IWR, *) 
WR I TE ( IWR, * ) 
WR I TE ( IWR, *) 
WRITE(IWR, *) 


'ASPECT RATIO=',AR 
'PLANFORM AREA (ft**2) = ',S 
'EFFICIENCY^', E 
'WEIGHT (lbs)-'.W 
'CDO= ', CDO 


V(ft/s> 
D( lb ) ' 


01 

01 

01 

01 

01 

01 

01 

01 

01 

01 


10 


C 

c 


WR I TE ( IWR, *) 

WRITE (IWR,*) ' 

CD 

1 = 1 

DO 10 V=10, 40 

X(I)=v 

CL=W/(. 001 19*(V**2. )*S) 

n D ^ C 59 +(CL ** 2)/<3 141 59*E*AR ) 
D=0. 001 19*CD*V**2. *S 
PR I=D*V 

Y(I, 1 )=PRI*1. 356 
PR=PR 1*1 356 

CONTfNii WR ’*’ V ' PR 'CL-CD,D 

GO TO 72 


PR ( W ) 


CL 


C 

C 


C 

C 


I0PT=-01 1 
ND=31 
NF= 1 

CAL}- Xr L0T( I0PT ' x ’ Y, 31, ND, NF) 

CALL TLABELC 'RELATIVE VELOCITY (ft/s)', 
CALL ^ITLE( 'PRELIMINARY POWER REQUIRED 


^POW|R^ REQUIRED (watts) 


- This 


WRITE (IWR,*) 

statement will never be executed 

IWR ' * ) ' V(ft/s) 

. EAL ( W ) pq ' 

DO 98 V=10, 40 
DO 99 N=4700, 9000, 500 
AN=N/60 
AJ=V/( AN*DI ) 

Pr- + 62 7 I?X +2 ; 57 64*AU-2. 7452*AJ**2. + 
pp"Iq A 1418*N-1. 1 731E-5*N**2. 

PA=PG?EP 53+ O0739 * N ~ 8 5894E-6*N**2. 
PAREAL=PA*1. 356 

contInue 1 ' * >v ' N ' EP ’ PAREAL ' p e 
CONTINUE 
CONTINUE 
CONTINUE 
STOP 
END 

Label 77 has not been referenced 


RPM 


EP 


01 

02 

02 

02 

02 

02 

02 

02 

02 

02 

01 


86713*AJ**3. 


99 

98 

77 

72 


END OF COMPILATION CLOCKED 


. 448 SECONDS 
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Appendix B-l 

Wing Analysis and TKSoIver! Code 


The wing is assumed to structurally consist of only a main spar as configured in the figures 
below. 


CM> V’EW 



T 1 * ii 3 
x * 12 

Io(r/ + Ad 1 ) 



* L/2 


For this analysis the maximum flight load factor is used to determine the loading condition. 
This is n=1.5. Lift force then equals n * W, where W is the estimated weight of 2.65 lbs. 
This force is distributed equally on both wings with a resultant at the midpoint of each side 
of the wing. This force creates a moment at the root of the wing given by 
M x = F * d = L/2 * b/4. The other force acting on the wing is due to drag. Drag force is 
neglected because (1) it is small compared to the lift force and (2) the moment of inertia 
about a y axis for the entire wing would actually take into account the effects of the leading 
and trailing edges. These edges, spread 8.5” apart greatly increase the Iy of the wing 
giving a small axial stress due to the drag force. The effect of the leading and trailing edges 
are neglected in analysis of the axial stress due to lift force because the contribution to I x of 
the leading and trailing edges is insignificant compared to that of the main spar. After 
moments have been calculated at the root the moment of inertia is calculated knowing that 
I x ' for a rectangle = 1/12 b * h3 (Reference 1). Using parallel axis theorem 
(lx = lx’ + Ad^, also Reference 1) the moment of inertia can be calculated for the entire 
spar. The above analysis has been coded and placed into TKSoIver as can be seen on the 
following page. Weight of the wing is calculated as well in this code, by accepting inputs 
of the densities of the various components as well as their cross sectional areas. Weight is 



then calculated by multiplying the length of the spar, leading edge, etc by their cross 
sectional area to get a volume and then multiplying by density to get a mass. 



S Rule 


* lp-(1/12*bp*t A 3)*2 

* fl-L/2*(b/1) 

* t-c*tc 

* flc-bc*hc 

* lo-I/12*(bc*hc A 3) 

* d-t/2-hc/2 

f lx»(lo+flc*d A 2)*2+lp 

* y-t/2 

* Slg-(n*y)/lx 

* Slgfa- Slg*l ,2 

* Qua I - Slg-Sigiax 

* U8par*flc*b*2*Rho$par 

* Q1«flc*d 

* iu1-(L*Qt)/(lx*t»eb) 

* H2«((t/2-hc)*t«b) 

* fltot-R2+flc 

* d2«(flc*d*R2*((t/2-hc)/2))/fltot 

* Q2«fltot*d2 

* Tau2-(L*Q2)/(lx*tteb) 

* Srlb-(t*c)/2 

* Nu«rib-b/4 

* Urlb-Srlb*trlb*fihorlb*Nu«rib 

* Ule-Sle*b*Rhole 

* Ute«Ste*b*Rhote 

* Ucover-b*c*2*fihocover 

* TotU»-Udpar+Ur I b*U I 8+Ut e+Ucover+Upane 1 8 

* S«lng«b*c 

* RR-b A 2/S*lng 

* Upane 1 3-bp*t *b*2*nhopane I 

* FS*Sig«ax/Sig 


Appendix B-2 

Fuselage Analysis 


The fuselage is configured as seen below. 




>> 



My = 9.1 lb Mx = 10.4 

y = 1 in x = 0*5 

Iy = 0.0365 in 4 lx * 0.1927 in 

The load on the fuselage is the load of the horizontal and vertical ails » laced “ 

end of the fuselage. The load carrying beams are assumed to be cann «« ^ ftom ^ 

wing box. The fuselage supports 0.35 lbs from fe > . 10 . 4 i„ lb. Using 

horizontal tail. This gives bending moments of My -9.1 in lb and 

the same formulas cited in Appendix HI,. 0.1927 in* and Iy = 0.0365 in . 
there are two components to the load the stress calculanon becomes: 

Axial Stress = Mu »x + MzJ-X 

Iy 

This equation is found from Allen and Haisletis toatoa '9 Aerospace Structu ral 
Analysis text. This yields a resultant axial stress of 178.63 pst. 


HORIZONTAL TAIL: 


The horizontal tail is analyzed in much the same way. Lift force on the tail is 0.4 lbs, 
found in stability and control analysis. This force is split between the two sides of the 
wing and placed at the midpoint of each side. It will be assumed that the leading and 
trailing edge of the non-movable area of the tail carry all of the load. Internal structure is 
non-load carrying and meant only to maintain the wing shape. Analysis follows as below. 
The stress in the members is found to be 365.4 psi. 









* 


Si 0e vVeuj 


M x = 1.9 in lb I x = 0.00065 in 4 
A = 0.0625 in 2 y= 0.125 in 


ORIGINAL PAGE IS 
0F POOR QUALITY 



Appendix B-3 

Empennage Analysis 


VERTICAL TAIL: 

Again a simple cantilever beam analysis is used Hie maximum load on the vertical tail is 
0.35 lbs as determined by the stability and control group. The vertical tail is configured as 

follows: 






e > jo view 


M = 1.75 in lb 
A = 0.0625 in 2 



v'/eW 


I x = 0.0026 in 4 
y = 0.25 in 


It is assumed that the load acts halfway up the tail giving a bending moment of M - L * h/2 
= 1.75 in lb. This moment will transfer direedy in to the pins which can be analyzed as 
cantilevered beams. See above diagrams. Calculating as in Appendix B-l axial stress is 
found to be 168 psi. It must also be ensured that the main area of the tail will not bre 
away from the seat. Assume that the leading and trading edge wiU carry all of the load. By 
the below analysis axial stress is 336 psi. 


t 1 

— — “X 

/ 

* * Tof \jia> 



M x = 1.75 in lb Ix - 0.00065 in 4 
A =0.0625 in 2 y = 0.125 in 




ORIGINAL PAGE IS 
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Appendix B-4 

Landing Gear Computer Code and Data 


This computer code makes use of the cantilever beam analysis of the strut of the landrng 
gear nte code varies the loading on the beam, the material of the beam and the radius 
length of the beam. This simple code provided data, the corresponding deflection ‘<°r the 
particular load, radius and length, to be analyzed for the sizing of the gear. The code an 
samples of the data provided are presented here in the code. 



Ul u l*j *■■■•■ 


c 

c 

c 


program tpaoe 

i' j R I t TEN 3/ - ERRy ZUROVCHAK 
AE441 DESIGN JVCHAK 

PEAL P, DELTA . : P . • 

j ^ •* EN~£P ;- r£T i »orrr ' 


5 WR I TE ( l , ♦ ; 
R E A D ( 1 , -* ) l 


‘ c ; A 


r H ( I r, J ' 


READ? 1 1 *?E “ N ' " ,K 4 ^C'ULUS OF ELAST I CITY ( p51 ) •- 

7r, . " ' • p 3 h ,• ' s? d r 1 'i? °i EUstic i <?y = ' , E 

3 i . 5 , loniin ) ' 


- . =*f. O * A TAN 1 

r re* ' - ; 

v t ■ A n - > : 


cc ; 


*. In A , -< ,» R 3 


i - w’ O n ■ _ 

t/ 0 DO J = ■ j. , ; * ” - . 

1= ■ P I *R-*-*4 * / 4 

r'CST*— ~ - 

■■" P * < ■- C J. l\ fr . - <3 r J 


3 Q 


0 \; T in* - 


I F * £, 

-OR MAT:. • r, 
3 ~GP 


p c*! T I 


rr -* 


5H TQ CONTINUE, ENTER 1 
3 X , Fit 3) 


C M 


'OT ENTER 


ORIGINAL PAGE IS 
of poor quality 
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Appendix C 


The information found in this Appendix is the program code and output, used to find the 
such as taper, apect ratio, and wing twist. 
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UNIVERSITY FTN77 VER 
OPTIONS: 
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0020 
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0029 

0030 

0031 

0032 

0033 

0034 

0035. 01 

0036. 01 

0037. 01 

0038. 01 

0039. 01 

0040. 02 

0041. 02 

0042. 02 

0043. 02 

0044. 02 

0045. 02 

0046. 02 

0047. 02 

0048. 02 

0049. 02 

0050. 02 

0051. 02 

0052. 02 

0053. 02 

0054. 02 

0055. 03 

0056. 03 

0057. 03 

0058. 03 

0059. 03 

0060. 03 
006 1 . 03 
0062. 03 

0063. 03 

0064. 04 

0065. 04 

0066. 04 

0067. 03 

0068. 02 

0069. 02 

0070. 02 

0071. 02 

0072. 02 

0073. 02 

0074. 02 i 

0075. 02 
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<USER2>S474988736>STU0ENT>FMECH. F77 

MOM A n i if-/' 


program fmech 
13 FEBRUARY J 1989 FER L ' NEU,1ANN 


NOSILESV E wJg? I I g,^ S ti 0 ^« R NOOFFSET LG( 


THIS PROGRAM 
AIRFOIL. 


WILL FIND THE FOURIER COEFFICIENTS FOR 


AN 


«lVcl j a n ' r' q' a/ 

REAL CL?GAMfe(S)fvELrTAU A sToMiT2! 4 it?i C 2' CT > C <4 ) , MU<4) . LCS 
REAL GANG (4), ALFA (4), 1(4 f l^TUPTA/ff 1 ^', ? UM < CDi * RATI0<4), AZL 
REAL CLS(8),DEGrx, P 4' 4 >* t HETA<4),Z(4),ANS<4>, IT# IR 
INTEGER F C 4 > , Y(4) 


;ENTER VARIABLE VALUES IN DEGREES. 


'WHERE 

'ENTER 


l 0 FOP U T^ N ^J^fk. PR0 © RA M TO 
1 FOR THE SCREEN AND 6 FOR 


RUN?' 

THE PRINTER. 


'THE 


VELOC I TY" 4 , VEL, ' f t/s ' 
ANGLE OF ZERO LIFT* ', AZL» 


deg. 


'THE 

3 


LIFT CURVE SLOPE"', LCS, 


per deg 


'THE ASPECT RATIO" ',AR 


'THE 

'WHAT 


WRITEC1, *) 

WRITEC1, *) 

WRITE(1, *) 

WRITEd, *) 

READ(1, *) Q 
VEL=25. 0 
WRITECG, *) 

AZL=-2. 0 
WRITECG,*) 'THE 
AL0=AZL/57. 3 
LCS"0. 09 
WRITECG, *) 

CLA=LCS*57. 

AR" 1 2. 00 
WRITECG, *) 

S=6. 00 
WRITE CG, * ) 

WRITEd, *) 

READC1,*) ia 
DO 99 KK"1, 15 
IR=KK 
IT"KK 

ALFAR=IR/57. 3 
ALFAT=IT/57. 3 
DO 60 H" 1 , 1 
TR(H) = 1. 0 
B=SQRTCAR*S) 

WRITECG, * ) ' ' 

writISq' J J ******************** ' 

WRITECG,*) 'b*',B, ' f eet' 

CT"TR C H ) #CR 
WRITECG,*) 'C t" ', CT 

WRITEC1,*) ' ' 


^Vc^rui AREA»',S, 'square ^eet' 

S THE INCIDENCE ANGLE IN DEGREES?' 


DO 


10 J=1 , N 

R*J*22. 5 

JHETAC J)= R/ C57. 3) 
OANoTSlii/S&STJJf S? < T HETA( J) , > 




MUCJ), ALFA C J ) 


& J W.’S ALFA(J? 

L=2*I-1 

CONTINUE > ~ SIN(L * THETA( J) >* (1 - 0+MUCJ)*L/SIN C THETA C J) ) ) 
CONTINUE 

DEG=ALFA C 1 ) *57 3 

WRITECG,*) 'ANGLE OF ATTACK = ', KK, ' DEGREES' 


CALL INVRTCA, N, 4, F, Y) 

WRITEuSo ;-ANS.4,4,1.4.4.4> 


WRITECG,*) 'THE COEFFICIENTS ARE. 
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Appendix D 


Appendix D is the derivation for the roll moment created by a dihedral angle and rudder 
input First of all, the change in angle of attack that the wing sees as a result of the dihedral 
angle is derived. Then, this is used to develop the equation few the roll moment. 
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